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Abstract: We present a general, rigorous theory of partition function zeros for lat- 
tice spin models depending on one complex parameter. First, we formulate a set of 
natural assumptions which are verified for a large class of spin models in a compan- 
ion paper [5]. Under these assumptions, we derive equations whose solutions give the 
location of the zeros of the partition function with periodic boundary conditions, up to 
an error which we prove is (generically) exponentially small in the linear size of the 
system. For asymptotically large systems, the zeros concentrate on phase boundaries 
which are simple curves ending in multiple points. For models with an Ising-like plus- 
minus symmetry, we also establish a local version of the Lee-Yang Circle Theorem. 
This result allows us to control situations when in one region of the complex plane the 
zeros lie precisely on the unit circle, while in the complement of this region the zeros 
concentrate on less symmetric curves. 
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1. Introduction 

1.1. Motivation. One of the cornerstones of equilibrium statistical mechanics is the 
notion that macroscopic systems undergo phase transitions as the external parameters 
change. A mathematical description of phase transitions was given by Gibbs [17] who 
characterized a phase transition as a point of non-analyticity in thermodynamic func- 
tions, e.g., the pressure. This definition was originally somewhat puzzling since actual 
physical systems are finite, and therefore their thermodynamic functions are manifestly 
real-analytic. A solution to this contradiction came in two seminal papers by Yang and 
Lee [25,42], where it was argued that non-analyticities develop in physical quantities 
because, as the system passes to the thermodynamic limit, complex singularities of 
the pressure pinch the physical (i.e., real) domain of the system parameters. Since the 
pressure is proportional to the logarithm of the partition function, these singularities 
correspond exactly to the zeros of the partition function. 

In their second paper [25], Lee and Yang demonstrated the validity of their theory in 
a particular example of the Ising model in a complex magnetic field h. Using an induc- 
tion argument, they proved the celebrated Lee- Yang Circle Theorem which states that, 
in this model, the complex-e'' zeros of the partition function on any finite graph with 
free boundary conditions lie on the unit circle. The subject has been further pursued by a 
number of authors in the following fifty years. Generahzations of the Lee- Yang theorem 
have been developed [26, 31, 32, 35] and extensions to other complex parameters have 
been derived (for instance, the Fisher zeros [14] in the complex temperature plane and 
the zeros of the (/-state Potts model in the complex-g plane [40,41]). Numerous papers 
have appeared studying the partition function zeros using various techniques including 
computer simulations [10,20,22], approximate analyses [21,24,29] and exact solu- 
tions of ID and 2D lattice systems [8, 9, 12, 18, 27, 28, 38, 39]. However, in spite of this 
progress, it seems fair to say that much of the original Lee- Yang program — namely, 
to learn about the transitions in physical systems by studying the zeros of partition 
functions — had remained unfulfilled. 

In [2], we outlined a general program, based on Pirogov-Sinai theory [6, 33, 34,43], 
to determine the partition function zeros for a large class of lattice models depending 
on one complex parameter z. The present paper, and its companion [5], give the math- 
ematical details of that program. Our results apply to a host of systems with first-order 
phase transitions; among others, they can be applied to field-driven transitions in many 
low-temperature spin systems as well as temperature-driven transitions — for instance, 
the order-disorder transition in the q'-state Potts model with large q or the confinement 
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Higgs transition in lattice gauge theories. We consider lattice models with a finite num- 
ber of equilibrium states that satisfy several general assumptions (formulated in detail 
below). The validity of the assumptions follows whenever a model can be analyzed us- 
ing a convergent contour expansion based on Pirogov-Sinai theory, even in the complex 
domain. In the present work, we study only models with periodic boundary conditions, 
although — with some technically involved modifications — our techniques should allow 
us to treat also other boundary conditions. 

Under our general assumptions, we derive a set of model-specific equations; the 
solutions of these equations yield the locations of the partition function zeros, up to 
rigorously controlled errors which are typically exponentially small in the linear size of 
the system. It turns out that, as the system size tends to infinity, the partition function 
zeros concentrate on the union of a countable number of simple smooth curves in the 
complex ^-plane. Another outcome of our analysis is a local version of the Lee- Yang 
Circle Theorem. Whereas the global theorem says that, for models with the full Ising 
interaction, all partition function zeros lie on the unit circle, our local theorem says that 
if the model has an Ising-like symmetry in a restricted region of the complex z-plane, 
the corresponding portion of the zeros lies on a piece of the unit circle. In particular, 
there are natural examples (see the discussion of the Blume-Capel model in [2]) where 
only some of the partition function zeros lie on the unit circle, and others lie on less 
symmetric curves. Our proof indicates that it is just the Ising plus-minus symmetry 
(and a natural non-degeneracy condition) that makes the Lee- Yang theorem true, which 
is a fact not entirely apparent in the original derivations of this result. 

In addition to being of interest for the foundations of statistical mechanics, our re- 
sults can often be useful on a practical level — even when the parameters of the model 
are such that we cannot rigorously verify all of our assumptions. We have found that 
our equations seem to give accurate locations of finite-volume partition function zeros 
for system sizes well beyond what can be currently achieved using, e.g., computer- 
assisted evaluations of these partition functions (see [2] for the example of the three 
dimensional 25-state Potts model on 1000 sites). Our techniques are also capable of 
handling situations with more than one complex parameter in the system. However, the 
actual analysis of the manifolds of partition function zeros may be technically rather 
involved. Finally, we remark that, in one respect, our program falls short of the ultimate 
goal of the original Lee- Young program — namely, to describe the phase structure of 
any statistical-mechanical system directly on the basis of its partition function zeros. 
Instead, we show that both the location of the partition function zeros and the phase 
structure are consequences of an even more fundamental property: the ability to rep- 
resent the partition function as a sum of terms corresponding to different metastable 
phases. This representation is described in the next section. 



1.2. Basic ideas. Here we will discuss the main ideas of our program, its technical diffi- 
culties and our assumptions in more detail. We consider spin models on l/, with d>2, 
whose interaction depends on a complex parameter z. Our program is based on the fact 
that, for a large class of such models, the partition function Z'^^'^ in a box of side L and 
with periodic boundary conditions can be written as 



r 




(1.1) 



m=l 
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Here gi, . . . , are positive integers describing the degeneracies of the phases 1, . . . , r, 
the quantities /i , . . . , are smooth (but not in general analytic) complex functions of 
the parameter z which play the role of metastable free energies of the corresponding 
phases, and f{z) = mini<„j<r Re/„i(z). The real version of the formula (1.1) was 
instrumental for the theory of finite-size scaling near first-order phase transitions [7]; 
the original derivation goes back to [6]. 

It follows immediately from (1.1) that, asymptotically as L tends to infinity, Z^""^ = 
requires that Re/m(z) = Re/m(z) = f{z) for at least two distinct indices m and rh. 
(Indeed, otherwise the sum in (l.I) would be dominated by a single, non- vanishing 
term.) Therefore, asymptotically, all zeros of Z'^' concentrate on the set 

^ = {z: there exist m ^ m with Kefm{z) = Refff^{z) = f{z)}. (1.2) 

Our first concern is the topological structure of Let us call a point where Re/m(z) = 
f{z) for at least three different m a multiple point; the points z € that are not 
multiple points are called points of two-phase coexistence. Under suitable assumptions 
on the functions /i, . . . , /r, we show that is a countable union of non-intersecting 
simple smooth curves that begin and end at multiple points. Moreover, there are only a 
finite number of multiple points inside any compact subset of C. See Theorem 2.1 for 
details. 

The relative interior of each curve comprising consists entirely of the points of 
two-phase coexistence, i.e., we have Re/„j(z) = Re/m(z) = f{z) for exactly two 
indices m and to. In particular, the sum in (1.1) is dominated by two terms. Supposing 
for a moment that we can neglect all the remaining contributions, we would have 

Z^iz) = + g^,e-/-(^)^^ (1.3) 

and the zeros of Z'^^ would be determined by the equations 

Re/m(z) = R&fih(z) + L^'^\og{qmlqfh) 
Imfraiz) = Im,fe(z) + {2£ + l)nL-\ 

where f is an integer. The presence of additional terms of course makes the actual zeros 
only approximate solutions to (1.4); the main technical problem is to give a reasonable 
estimate of the distance between the solutions of (1.4) and the zeros of Z"^' . In a neigh- 
borhood of multiple points, the situation is even more complicated because there the 
equations (1.4) will not be even approximately correct. 

It turns out that the above heuristic argument carmot possibly be converted into a 
rigorous proof without making serious adjustments to the initial formula (1.1). This is 
a consequence of subtle analytic properties of the functions /,„. For typical physical 
systems, the metastable free energy is known to be analytic only in the interior of 
the region 

ym = {z:'R-Qfm{z) = f{z)]. (1.5) 

On the boundary of one expects — and in some cases proves [15, 19] — the existence 
of essential singularities. Thus (1.1) describes an approximation of an analytic function, 
the function Z"^' , by a sum of non-analytic functions, with singularities appearing pre- 
cisely in the region where we expect to find the zeros of Z^' ! It is easy to construct 
examples where an arbitrarily small non-analytic perturbation of a complex polynomial 
with a degenerate zero produces extraneous roots. This would not be an issue along the 
two-phase coexistence lines, where the roots of Z'^^^ turn out to be non-degenerate, but 
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we would not be able to say much about the roots near the multiple points. In short, we 
need an approximation that respects the analytic structure of our model. 

Fortunately, we do not need to look far to get the desirable analytic counterpart of 
(1.1). In fact, it suffices to modify slightly the derivation of the original formula. For 
the benefit of the reader, we will recall the main steps of this derivation: First we use 
a contour representation of the model — the class of models we consider is character- 
ized by the property of having such a contour reformulation — to rewrite the partition 
function as a sum over the collections of contours. Then we divide the configurations 
contributing to Z^'^ into r + 1 categories: Those in which all contours are of diameter 
smaller than, say, L/3 and in which the dominant phase is to, where m = 1, .... r, and 

those not falling into the preceding categories. Let Zm"^ be the partial partition function 
obtained by summing the contributions corresponding to the configurations in the TO-th 
category, see Fig. 1 . It turns out that the error term is stiU uniformly bounded as in ( 1 . 1), 
so we have 

Zl"{z) = ^(^) + 0(e— ^e-^'(^)^'), (1.6) 

but now the functions Zm\z) are analytic, and non-zero in a small neighborhood 
of S^m- (However, the size of the neighborhood shrinks with L — > oo, and one of 
the challenges of using the formula (1.6) is to cope with this restriction of analyticity.) 
Moreover, writing 

Z(f)(.) = g„e-/i^'W^^ (1.7) 

and using the contour representation, the functions fm^ can be expressed by means 
of convergent cluster expansions [11, 23]. In particular, they can be shown to converge 
quickly to the functions as L ^ oo. 

In this paper, we carry out the analysis of the partition function zeros starting from 
the representation (1.6). In particular, we formulate minimal conditions (see Assump- 
tions A and B in Sect. 2) on the functions fm^ and the error terms that allow us to 
analyze the roots of Z'^^ in great detail. The actual construction of the functions fm^ 
and the proof that they satisfy the required conditions is presented in [3,4] for the q- 
state Potts model with one complex external field and q sufficiently large, and in [5] for 
a general class lattice models with finite number of equilibrium states. 



1.3. Discussion of assumptions and results. Here we will describe our main assump- 
tions and indicate how they feed into the proofs of our main theorems. For consistency 
with the previous sections, we will keep using the functions and fm^ even though 
the assumptions will actually be stated in terms of the associated exponential variables 

= e-/-^^) and Ci^)(^) = e"/- (1.8) 

The first set of assumptions (Assumption A, see Sect. 2.1) concerns the infinite-volume 
quantities /„, and is important for the description of the set of coexistence points 
The functions are taken to be twice differentiable in the variables x = Re2 and y = 
\mz, and analytic in the interior of the set <5^m- If, in addition, f{z) — min^ ^sfm is 
uniformly bounded from above, good control of the two-phase coexistence curves is ob- 
tained by assuming that, for any distinct to and to, the difference of the first derivatives 
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+ ^ + + ^ 
+ + + + + 



Fig. 1. Schematic examples of configurations, along with their associated contours, which contribute to dif- 
ferent terms in the decomposition in ( 1 .6). Here we have a spin model with r = 3 equilibrium phases denoted 
by +, — and 0. The configuration on the left has all contours smaller than the cutoff — which we set to L/3 

where L is the side of the box — and will thus contribute to because + is the external phase for all 

external contours. The configuration on the right has long contours and will be assigned to the error term. 



of fm and fm is uniformly bounded from below on n Finally, in order to dis- 
cuss multiple coexistence points, we need an additional non-degeneracy assumption on 
the derivatives of the functions /,„ for the coexisting phases. Given these assumptions, 
we are able to give a very precise characterization of the topology of the coexistence 
set Sf, see Theorem 2.1. 

The second set of assumptions (Assumption B, see Sect. 2.2) is crucial for our re- 
sults on the partition function zeros, and is formulated in terms of the functions fm^. 
These will be taken to be analytic with a uniform upper bound on the first r derivatives 
in an order-(l/iy) neighborhood of the sets .ym- In this neighborhood, /if •* is also as- 
sumed to be exponentially close to /„, with a lower bound on the difference of the first 
derivatives for any pair fm^ and in the intersection of the corresponding order- 
(1/L) neighborhoods. Finally, we need a bound on the error term and its derivatives in 
an approximation of the form (1 .6) where the sum runs only over the dominating terms, 
i.e., those m for which z lies in the order-(l / L) neighborhood of J^^. 

Combining Assumptions A and B, we are able to prove several statements on the 
location of the partition function zeros. We will start by covering the set of available 
2;-values by sets with a given number of stable (or "almost stable") phases. The cov- 
ering involves three scale functions, cjl, 7l and pL which give rise to three classes of 
sets: the region where one phase is decisively dominating the others (more precisely, 
the complement of an L^'^o^l -neighborhood of the set ^^), a 7^ -neighborhood of sets 
with two stable phases, excluding a 71, -neighborhood of multiple points, and the pL- 
neighborhoods of multiple points. As is shown in Proposition 2.6, for a suitable choice 
of sequences ojl, 7^, and pL, these three sets cover all possibilities. 

In each part of the cover, we will control the zeros by a different method. The re- 
sults of our analysis can be summarized as follows: First, there are no zeros of 
outside an L~''cji-neighborhood of the set f^. This claim, together with a statement 
on the maximal possible degeneracy of zeros, is the content of Theorem 2.2. The next 
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theorem, Theorem 2.3, states that in a 7^ -neighborhood of the two-phase coexistence 
points, excluding a neighborhood of multiple points, the zeros of Z'^' are exponen- 
tially close to the solutions of (1.4). In particular, this implies that the zeros are spaced 
in intervals of order-L^'^ along the two-phase coexistence curves with the asymptotic 
density expressed in terms of the difference of the derivatives of the corresponding free 
energies — a result known in a special case already to Yang and Lee [42]; see Proposi- 
tion 2.4. The control of the zeros in the vicinity of multiple points is more difficult and 
the results are less detailed. Specifically, in the -neighborhood of a multiple point 
with q coexisting phases, the zeros of Z^' are shown to be located within a L"'*-'^/? 
neighborhood of the solutions of an explicitly specified equation. 

We finish our discussion with a remark concerning the positions of zeros of complex 
functions of the form: 

r 

Zn{z) = J2 (^m{z)Cm{z)'', (1.9) 

m—1 

where ai, . . . , and (i, . . . ,(r are analytic functions of z. Here there is a general 
theorem, due to Beraha, Kahane and Weiss [1] (generalized recently by Sokal |41]), 
that the set of zeros of Zn asymptotically concentrates on the set of z such that ei- 
ther am{z) = and |Cm(2:)| = maxfe \Ck{z)\ for some m = 1, ... ,r or \Cm{z)\ = 
\Cn{z)\ = maxfe |Cfe(2)| for two distinct indices m and n. The present paper provides a 
substantial extension of this result to situations when analyticity of Cm {z) can be guar- 
anteed only in a shrinking neighborhood of the sets where m is the "dominant" index. 
In addition, we also provide detailed control of the rate of convergence. 



2. Main results 

2.1. Complex phase diagram. We begin by abstracting the assumptions on the meta- 
stable free energies of the contour model and showing what kind of complex phase 
diagram they can yield. Throughout the paper, we will assume that a domain (Z C 
and a positive integer r are given, and use TZ to denote the set 7?. = {1, . . . , r}. For 
each z G ^,we let x = Rez and y = Imz and define, as usual, 

9^ = U£-^fy) d,= Ui,+ify)- (2.1) 

Assumption A. There exists a constant a > and, for each m & TZ, a function 
(m- G ^ such that the following conditions are satisfied: 

(1) The quantity C(^) = max^gK I Cm (2)! is uniformly positive in G, i.e., we have 

inf^e£?C(2) > 0. 

(2) Each function viewed as a function of two real variables x = Re 2; and y = 
Imz, is twice continuously differentiable on and it satisfies the Cauchy-Riemann 
equations ^fCml-^) = for all 2 e J^m. where 

ym = {z&G:\Cm{z)\=C{z)}- (2.2) 

In particular. Cm is analytic on the interior of .S^m- 

(3) For any pair of distinct indices m.n ^TZ and any z G S^m H S^n we have 



dzCrn{z) dzCn{z) 



C,m{z) Cn(^) 



> a. (2.3) 
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(4) If Q c 7e is such that | Q| > 3, then for any 2; e flmea 

^m(-2) = /- ( ^ ^ me Q, (2.4) 

are the vertices of a strictly convex polygon in C ~ M^. 

Remark 1. In (1), we assumed uniform positivity in order to simpUfy some of our later 

arguments. However, uniformity in can easily be replaced by uniformity on compact 
sets. Note that Assumptions A3-4 are invariant with respect to conformal transforma- 
tions of ^ because the functions involved in (2.3) and (2.4) satisfy the Cauchy-Riemann 
conditions. Also note that, by Assumption A3, the length of each side of the polygon 
from Assumption A4 is at least q ; cf Fig. 3. 

The indices m eTZ will be often referred to as phases. We call a phase m stable at z 
a z G S^rn, i-c, if \Cm{z) \ = C{z). For each z G ^ we define 

Q{z) = {men:\Cm{z)\=az)} (2.5) 

to be the set of phases stable at z.lfm,n G Q{z), then we say that the phases m and n 
coexist at z. The phase diagram is determined by the set of coexistence points: 

?f= y ^(m,n) with ^(m,n) = ^„ n =5^„. (2.6) 

If |Cm(^)| = C{z) for at least three distinct m G TZ, we call such z G ^ a multiple 
point. 

In the following, the phrase simple arc denotes the image of (0, 1) under a continu- 
ous and injective map while simple closed curve denotes a corresponding image of the 
unit circle {z G C: \z\ = 1}.A curve will be called smooth if it can be parametrized 
using twice continuously differentiable functions. 

Our main result concerning the topology of ^ is then as follows. 

Theorem 2.1. Suppose that Assumption A holds and let G ff be a compact set. 
Then there exists a finite set of open discs Di, D2, • • • , C ^ covering Si, such that 
for each k = 1, . . . ,i, the set s^/k = ^ n satisfies exactly one of the following 
properties: 

(1) -s^k = 0. 

(2) £/k is a smooth simple arc with both endpoints on dOk- Exactly two distinct 
phases coexist along the arc constituting s^k- 

(3) s^k contains a single multiple point Zk with Sk = \Q{zk)\ > 3 coexisting phases, 
and M; \ {zk} is a collection of Sk smooth, non-intersecting, simple arcs con- 
necting Zk to 9Dfe. Each pair of distinct curves from s^k \ {zk} intersects at a 
positive angle at Zk. Exactly two distinct phases coexist along each component 
ofs^k \ {zk}- 

In particular, = U^ee "^^^^^ ^ is a finite or countably -infinite collection of 
smooth simple closed curves and simple arcs which intersect each other only at the 
endpoints. 

Theorem 2.1 is proved in Sect. 3.2. Further discussion is provided in Sect. 2.4. 
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Fig. 2. An illustration of the sets ^iCQ) in the vicinity of a multiple point. The thick lines indicate the 
visible portion of the set of coexistence points Sf. Three phases, here labeled 1, 2 and 3, are stable at the 
multiple point. In (a), the three shaded domains represent the sets '^^({l}), '^e({2}) and "^^({3}), with the 
label indicated by the number in the box. Similarly, in (b) the three regions represent the sets '^e({l, 2}), 
'^e({2, 3}) and '^e({l, 3}). Finally, (c) contains only one shaded region, representing the set'We{{l, 2, 3}). 
The various regions '^e (Q) generously overlap so that their union covers the entire box. 



2.2. Partition function zeros. Next we will discuss our assumptions and results con- 
cerning the zeros of the partition function. We assume that the functions Z^' : ^ — > C, 
playing the role of the partition function in a box of side L with periodic boundary con- 
ditions, are defined for each integer L, or, more generally, for any L Gh, where L C N 
is a fixed infinite set. Given any m gTZ and e > 0, we use ^e(m) to denote the region 
where the phase m is "almost stable," 

^e{m) = {ze&: \Cm{z)\ > e-%{z)}. (2.7) 

For any Q C 7?., we also introduce the region where all phases from Q are "almost 
stable" while the remaining ones are not, 

^e(Q) = f] J'e{m)\\J ^e/2{n), (2.8) 

with the bar denoting the set closure. Notice that the function (m is non-vanishing 
on ye{m) and that Uqck '^e{Q) = ^, see Fig. 2. Note also that = 0, so we 

may assume that Q ^ for the rest of this paper. 

Assumption B. There exist constants k,t G (0, oo) and, for each m G TZ, a positive 

integer and a function (^m^ : y^/L {fn) ^ C such that for any L e L the following 
is true: 

(1) The function Z^" is analytic in 

(2) Each ' is non-vanishing and analytic in Furthermore, 



log 



< e"^^ (2.9) 
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and 



dzlog 



Cm{z) 



92 log 



C^\z) 

Cm{z) 



< e 



-tL 



(2.10) 



for all TO e 7?. and all z G c^k/l(to). (Here "log" denotes the principal branch of 
the complex logarithm.) 
(3) There exist constants d > 0, M < oo and Lq < oo such that for any L > Lq we 
have 

diC^\z) 



C.li\z) 



< M, 



whenever m € TZ,i= 1, . . . , r, and z G =5^k/l(to). In addition, 

d.(:^^\z) d,ci''\z) 



> a 



(2.11) 



(2.12) 



whenever TO, n e 7?. are distinct and z e fl ^K,/L{n). 

(4) There exist constants < oo, = 0, 1, . . . , r + 1, such that for any Q C 7?., the 
difference 

^qA^) = zt\z) - q^[c^,;:\z)f n.n) 

satisfies the bound 



meTl 



-tL 



(2.14) 



for all ^ = 0, 1, . . . , r + 1, uniformly in ^ e ^^^/z, (Q). 



Remark 2. In applications, will represent the degeneracy of the phase rn; thus we 
have taken it to be a positive integer. However, our arguments would go through even 
if we assumed only that all gm's are real and positive. It is also worth noting that 
in many physical models the partition function is not directly of the form required 
by Assumption B; but it can be brought into this form by extracting a multipUcative 
"fudge" factor F{z)^ , where F{z) ^ in the region of interest. For instance, in the 
Ising model with z related to the complex external field ft, by z = e'* we will have 



to take F{z) = z 
of 2 = 0. 



-1/2 



to make the partition function analytic in the neighborhood 



Our first theorem in this section states that the zeros of Z^'^{z) are concentrated 
in a narrow strip along the phase boundaries. In addition, their maximal degeneracy 
near the multiple points of the phase diagram can be evaluated. In accord with the 
standard terminology, we will call a point zq a k-times degenerate root of an analytic 
function h,{z) if h,(z) — g{z){z — z^)'' for some g{z) that is finite and non-zero in a 
neighborhood of zq. Recalling the definition (2.8) of the set "^(4(2), we introduce the 
shorthand 

^e=\J (^e/2(n)n^,/2(TO)) |J 'We{{m}) . (2.15) 



meTl 



An easy way to check the second equality in (2.15) is by noting that & \ '^^{{m}) can 
be written as the union Un nT^m =^e/2(^)- Then we have the following result. 
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Theorem 2.2. Suppose that Assumptions Al-3 and B hold and let k > be as in As- 
sumption B. Let {wl) bea sequence of positive numbers such that lol ^ cg. Then there 
exists a constant Lq < oo such that for L> Lq all roots ofZ'^' lie in '^l—'ul 
at most \TZ\ — 1 times degenerate. For each Q C TZ, the roots of in '^k/l{Q) <^^^ 
at most I Q| — 1 times degenerate. 

In other words, as L ^ oo, the zeros of Z'^^ asymptotically concentrate on the set of 
coexistence points $f . Notice that we exphcitly do not require Assumption A4 to hold; 
see Sect. 2.4 for further discussion. Theorem 2.2 is proved in Sect. 4.1. 

Our next theorem deals with the zeros of Z'^^'^ in the regions where at most two 
phases from TZ are "almost stable." It turns out that we have a much better control on 
the location of zeros in regions that are sufficiently far from multiple points. To quantify 
the meaning of "sufficiently far," we let 7l be a sequence of positive numbers (to be 
specified below) and, for any Q c 7?. with | Q| = 2 and any L > 0, let : -^^C^^ (Q) 
(0, oo) be a function defined by 

. . ^ if^e"^7J2)n^2«/L(Q), ... 

^'^'^ = \L''e-i-^^ otherwise. ^^''^^ 

(Clearly, 5l{z) depends on the index set Q. However, this set will always be clear from 
the context and so we will not make it notationally explicit.) Finally, given e > and 

2: G let ^c{z) denote the open disc of radius e centered at z. 

The exact control of the roots in two-phase regions is then as follows. 

Theorem 2.3. Suppose that Assumptions A and B hold and let i?^ be the set of all zeros 
of the function Z^'^{z) in G, including multiplicity. If m,n € TZ are distinct indices, 
let Q = {m, n}, and let QhiQ) be the set of the solutions of the system of equations 

qT^CA^)\-^in"''\Uz)l (2.17) 

i''Arg(C™(z)/C„(z)) =7rmod27r. (2.18) 

Let be such that 

liminf— ^>4d and limsupL'^-V < 2t, (2.19) 

L^oo logX i^oo 

and let Sl '■ ^71, (Q) (0, 00) be as defined in (2.16). Then there exist finite positive 
constants B, C, D, and Lq such that for any Q C TZ with \Q\ = 2 and any L > Lq 
we have: 

(1) For all z e n '^7i,(Q) with Iiij:,i^-d{z) C G, the disc Il)uj^-d{z) contains at 
least one root from Q^. 

(2) For all z G (1 '^^^{Q) with I}c5Liz){z) C G, the disc 3c5^{z){z) contains 
exactly one point from f2i{Q). 

(3) Forallz € ^2^(2)0'^-^^ (Q) withBcsL{z){z) C the disc'S3c5L(z){z) contains 
exactly one root from fi*i^. 

(4) Any two distinct roots of Z^^ in the set {2 e (Q) : "QQi^-diyZ) C ^} are at 
least BL^"^ apart. 
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Note that the first limit in (2.19) ensures that L'^5l{z) — > as L ^ oo throughout 
'^"iL (2) (for any Q C TZ with | Q| = 2). Thus 6l {z) is much smaller than the distance 
of the "neighboring" roots of (2.17-2.18). Theorem 2.3 is proved in Sect. 4.2. 

Theorem 2.3 allows us to describe the asymptotic density of the roots of Z'^" along 
the arcs of the complex phase diagram. Let m,n G TZhe distinct and let ^^(m, n) be as 

in (2.6). For each e > and each z € $^(to, n), let pm^'n (z) be defined by 

&iz) = ^\^2tr^Mz)\, (2.20) 

where fl ©^(-z)! is the number of roots of Z^j^^ in M^{z) including multiplicity. 
Since ^^(m, n) is a union of simple arcs and closed curves, and since the roots of (2.17- 

2. 18) are spaced within 0{L~'^) from each other, pm',n {z) has the natural interpretation 
of the approximate line density of zeros of Z"^^ along $f (m, n) . As can be expected from 

Theorem 2.3, the approximate density pm',n {z) tends to an explicitly computable Umit. 

Proposition 2.4. Let m,n e TZbe distinct and let pm^'n (z) be as in ( 2.20). Then the 
limit 

Pm,n{z) = lim lim {z) (2.21) 

ej.0 I/— >cx) 

exists for all z G $f(m, n) such that \Q{z)\ = % and 



Pm,n{z) — 2jj. 



dzCm{z) dzCn{z) 



Cm{z) Cn{z) 



(2.22) 



Remark 3. Note that, on the basis of Assumption A3, we have that Pm,n{z) > ct/ (27r). 
In particular, the density of zeros is always positive. This is directly related to the fact 
that all points z will exhibit a first-order phase transition (defined in an appropriate 
sense, once Imz ^ or Rez < 0) — hence the title of the paper The observation that the 
(positive) density of zeros and the order of the transition are closely related goes back 
to [42]. 

In order to complete the description of the roots of ZY', we also need to cover 
the regions with more than two "almost stable" phases. This is done in the following 
theorem. 

Theorem 2.5. Suppose that Assumptions A and B are satisfied. Let Zu be a multiple 
point and let Q— Q{zm) with q — \Q\ > 3. For each m G Q, let 

cPm{L) = L'^ AigCmizu) {mod2n) and = (2.23) 

(,m[ZM) 

Consider the set Ql (Q) of all solutions of the equation 

^{l)+l'^{z-zm)v^ _ (2.24) 



including multiplicity, and let {pl) bea sequence of positive numbers such that 

lim L'^pL = 00 but lim = 0. (2.25) 

L^oo L^oo 
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Define p'^ = Pl -\- L'^'^f^+^Z''). Then there exists a constant Lq < co and, for any L > 
Lq, an open, connected and simply connected set satisfying Dp^ (^m) C ^ C 
{zm) such that the zeros in f2f\^ are in one-to-one correspondence with the solu- 
tions in i7(Q)n^ andthe corresponding points are not farther apart than L~'*(i+i/9). 

Theorem 2.5 is proved in Sect. 4.4. Sect. 2.4 contains a discussion of the role of 
Assumption A4 in this theorem; some information will also be provided concerning the 
actual form of the solutions of (2.24). 

To finish the exposition of our results, we wiU need to show that the results of The- 
orems 2.2, 2.3 and 2.5 can be patched together to provide complete control of the roots 
of Z"^^, at least in any compact subset of ^. This is done in the following claim, the 
proof of which essentially reUes only on Assumption A and compactness arguments: 

Proposition 2.6. Suppose that Assumption A holds and let uj^, 7^ and pL be sequences 
of positive numbers such that u>l < "JlL'^, 7l 0, and pL 0. For each compact 
set Si d 0, there exist constants x = x(^) > and Lq = Lo{S) < 00 such that, if 
Pl > XI L, we have 

^i-.,,n^c y %AQ)^ U (2.26) 

Qcn zue^ 

|Q|=2 |Q(zm)|>3 

for any L > Lq. 

Note that in (2.26) we consider only that portion of 3i in '^l-'^wl^ since by The- 
orem 2.2 the roots of Z'^^^ are contained in this set. Note also that the conditions we 
impose on the sequences ul, 7l and pL in Theorems 2.1, 2.3 and 2.5 and Proposi- 
tion 2.6 are not very restrictive. In particular, it is very easy to verify the existence 
of these sequences. (For example, one can take both 7^ and p^ to be proportional to 
i"'' log i with suitable prefactors and then let = L'^jl-) 

2.3. Local Lee-Yang theorem. As our last result, we state a generalized version of the 
classic Lee- Yang Circle Theorem [25], the proof of which is based entirely on the exact 
symmetries of the model. 

Theorem 2.7. Suppose that Assumptions A and B hold. Let -\- and — be two selected 
indices from TZ and let ^ be an open set with compact closure Si d such that 
^r\{z: \z\ = 1} ^ ^. Assume that S is invariant under circle inversion z l/z*,and 

(1) zi\z) = zi\\iz*r, 

(2) C+{z) = C-{l/z*randq+ = q- 

hold for all z ^ S) and all L S L. Then there exists a constant Lq such that the following 
holds for all L > Lq: If the intersection of S with the set of coexistence points ^ is 
connected and if-\- and — are the only stable phases in S, then all zeros in S lie on the 
unit circle, and the number of zeros on any segment of Sir\{z: \z\ = 1} is proportional 

to L'^ as L 00. 

Condition (2) is the rigorous formulation of the statement that the -|- and — phases 
are related by 2 1/z* (or h <-> —h, when z = e^) symmetry. Condition (1) then 
stipulates that this synometry is actually respected by the remaining phases and, in par- 
ticular, by Z^"' itself. 
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Remark 4. As discussed in Remark 2, in order to satisfy Assumption B it may be nec- 
essary to extract a multiplicative "fudge" factor from the partition function, perform 
the analysis of partition function zeros in various restricted regions in C and patch the 
results appropriately. A similar manipulation may be required in order to apply Theo- 
rem 2.7. 

Here are the main steps of the proof of Theorem 2.7: First we show that the phase 
diagram in ^ falls exactly on the unit circle, i.e., 

= {z e ^: |z| = 1}. (2.27) 

This fact is essentially an immediate consequence of the symmetry between "-h" and 
"— ." A priori one would then expect that the zeros are close to, but not necessarily on, 
the unit circle. However, the symmetry of combined with the fact that distinct ze- 
ros are at least BL~'^ apart is not compatible with the existence of zeros away from the 
unit circle. Indeed, if ^; is a root of Z^'^, it is bound to be within a distance 0{e~'^^) of 
the unit circle. If, in addition, \z\ ^ 1, then the z ^ \/ z* symmetry implies that 1/z* 
is also a root of Z'^^\ again within 0(e^^^) of the unit circle. But then the distance be- 
tween z and 1/2:* is of the order e~^^ which is forbidden by claim (4) of Theorem 2.3 . 

This argument is made precise in the following proof. 

Proof of Theorem 2.7. We start with the proof of (2.27). Let us suppose that S! C and 
Q{z) C {+, — } for all z £ Invoking the continuity of C,± and condition (2) above, 
we have Q{z) = { + , -} for all z e ^r\{z: |z| = 1} and thus i?n {z: |z| = 1} C ^. 
Assume now that ^ n ^ \ {z : |z| = 1} 9^ 0. By the fact that D ^ is connected and 
the assumption that n{z: \z\ = 1} ^ 0, we can find a path zt € ^ f) e [-1,1], 
such that zt e ^n{z: |z| = 1} if i < and zt € ^ n ^ \ {z : |z| = 1} if t > 0. Since 
Q(zo) = { + ,—}, we know that there is a disc Dg(zo) C that contains no multiple 
points. Applying Theorem 2. 1 to this disc, we conclude that there is an open disc D with 
zo € ID C De(zo), such that 5f fl D is a simple curve which ends at i9ID). However, using 
condition (2) above, we note that as with Zt, also the curve t l/zj lies in ^ n ^, 
contradicting the fact that !^ n B is a simple curve. This completes the proof of (2.27). 

Next, we will show that for any zo & S' (1 {z: \z\ = 1}, and any S > 0, there 
exists an open disc ©^(zo) C such that the set fl ©^(-20) is a smooth curve with 
the property that for any z G ID'e(zo) with |z| 7^ 1, the line connecting z and 1/z* 
intersects the curve n 3e{zo) exactly once, and at an angle that lies between tt/2 — S 
and tt/2 + S. If zq lies in the interior of S^, this statement (with (5 = 0) follows trivially 
from (2.27). If zq is a boundary point of ^, we first choose a sufficiently small disc 
]D> 9 zo so that D c 1^ and, for all points in D, only the phases -|- and — are stable. Then 
we use Theorem 2.3 and (2.27) to infer that e can be chosen small enough to guarantee 
the above statement about intersection angles. 

Furthermore, we claim that given zq G ^ (1 {z: \z\ = 1} and e > such that 
ID)3£(zo) C ^ and Q(z) C {+, — } for all z € P3e(zo), one can choose L sufficiently 
large so that 

»2e(2o) n ^^i-.^, c ({+, -}) n ^2./Li{+, -})• (2.28) 

To prove this, let us first note that, for 7^ < 2k/L, the right hand side can be rewritten as 



(2.29) 
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Next, by the compactness of D2e(zo) and the fact that no m G 7?. different from ± is 
stable anywhere in 'S]>^^{zq), we can choose Lq so large that ^ ^2e{zo) = 

for all L > Lq and all m ±. Using the closure of D2£(2;o) in place of the set ^ in 
(2.26), we get (2.28). 

We are now ready to prove that for any zq e ^n{z: \z\ = 1}, there exist constants 
e > and Lq such that all roots of ^£,"in (zo)n^ lie on the unit circle. To this end, let 
us first assume that e has been chosen small enough to guarantee that (1 — e)^^ < l + 2e, 
^3e{zo) C ff, Q{z) C {+, -} for all e B3e{zo), and ^ n na^izo) is a smooth curve 
with the above property about the intersections angles, with, say, 5 = 7r/4. Assume 
further that L is chosen so that (2.28) holds and e > max{C5L{zo), BL~'^), where C 
and B are the constants from Theorem 2.3. 

Let z e ^e{zo) n ^ be a root of If L is so large that Theorem 2.2 applies, we 
have z G '^^l-'^ul and thus 5l{z) = e~'^^ in view of (2.28). By Theorem 2.3, there 
exists a solution z to (2.17-2.18) that lies in a Ci5L(z)-neighborhood of z, implying 
that z has distance less than C6l{z) from D2<;(zo) H (Here we need that g+ = gr_ 
to conclude that z G '^.) Suppose now that \z\ ^ 1. Then the condition (I) above 
implies that z' = {z*)~^ is a distinct root of Z^'^ in ^. Moreover, if e is so small 
that (1 - e)"i < 1 + 2e, then z' € ^l~<^^l ^ "^^eXzo) and bhiz') also equals e"'^'^, 
implying that z' has distance less than C6h{z) from 035(20) H ^. Since both 2 and z' 
have distance less than C5l{z) fromD3e(2;o) fl^^, and the curve Us,; (2:0) fl^ intersects 
the line through z and z' in an angle that is near tt/2, we conclude that \z — z'\ < 
2\f2Ce~'^^ which for L sufficiently large contradicts the last claim of Theorem 2.3. 
Hence, z must have been on the unit circle after all. 

The rest of the argument is based on compactness. The set ^ n {z: |z| = 1} is 
compact, and can thus be covered by a finite number of such discs. Picking one such 
cover, let be the complement of these disc in '3. Then the set Si' is a finite distance 
away from ^ and thus 2l' fl '^l-'^wl = for L sufficiently large. From here it follows 
that for some finite Lq < 00 (which has to exceed the maximum of the corresponding 
quantity for the discs that constitute the covering of ^ fl {2 : \z\ = 1}), all roots of Z^'^ 
in he on the unit circle. □ 



2.4. Discussion. We finish with a brief discussion of the results stated in the previ- 
ous three sections. We will also mention the role of (and possible exceptions to) our 
assumptions, as well as extensions to more general situations. 

We begin with the results on the complex phase diagram. Theorem 2. 1 describes the 
situation in the generic cases when Assumptions A1-A4 hold. We note that Assump- 
tion A3 is crucial for the fact that the set ^ is a collection of curves. A consequence 
of this is also that the zeros of Z^''' asymptotically concentrate on curves — exceptions 
to this "rule" are known, see, e.g., [36]. Assumption A4 prevents the phase coexistence 
curves from merging in a tangential fashion and, as a result of that, guarantees that 
multiple points do not proliferate throughout ff. Unfortunately, in several models of in- 
terest (e.g., the Potts and Blume-Capel model) Assumption A4 happens to be violated 
at some z for one or two "critical" values of the model parameters. In such cases, the 
region i? has to be restricted to the complement of some neighborhood of z and, inside 
the neighborhood, the claim has to be verified using a refined and often model-specific 
analysis. (It often suffices to show that the phase coexistence curves meeting at z have 
different curvatures, which amounts to a statement about the second derivatives of the 
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Fig. 3. An illustration of the situation around a quadruple point. Here , . . . , t)| are the complex conjugates 
of the quantities from (2.4) and qi = 92 = 93 < 94 ■ (The quadruple point lies at the common tail point of the 
vectors v^, . . . ,vl.) The dashed lines indicate the asymptotes of the "strings" of zeros sufficiently far — on 
the scale L~'* — from the quadruple point. Note the lateral shift of these lines due to the fact that 94 > qi , 93. 
The picture seems to suggest that, on the scale L~'^, the quadruple point spUts into two triple points. 



functions log (m{z).) Examples of such analysis have appeared in [2] for the Blume- 
Capel model and in [4] for the Potts model in a complex external field. 

Next we will look at the results of Theorems 2.2 and 2.3. The fact that the roots 
of Z^'^ are only finitely degenerate is again independent of Assumption A4. (This is of 
some relevance in view of the aforementioned exceptions to this assumption.) The fact 
that, in the cases when all g^'s are the same, the zeros shift only by an exponentially 
small amount away from the two-phase coexistence lines is a direct consequence of our 
choice of the boundary conditions. Indeed, the factor e~^^ in (2. 1 6) can be traced to the 
similar factors in (2.9) and (2.14). For strong (e.g., fixed-spin) boundary conditions, we 
expect the corresponding terms in (2.9) and (2.14) to be replaced by 1 /L. In particular, 
in these cases, the lateral shift of the partition function zeros away from the phase- 
coexistence lines should be of the order 1/L. See [44] for some results on this problem. 

Finally, let us examine the situation around multiple points in some detail. Theo- 
rem 2.5 can be given the following geometrical interpretation: Let Zm be a multiple 
point. Introducing the parametrization ^ ~ (z ~ zm)L'^, we effectively zoom in on the 
scale L^"*, where the zeros of Z^'^ are well approximated by the roots of the linearized 
problem (2.24) with Q = Q{zm)- Let us plot the complex conjugates of the loga- 
rithmic derivatives Vm (see (2.23)), me Q, as vectors in M?. By Assumption A4, the 
vectors are the endpoints of a convex set in C ~ K^. Let vl, . . . ,v*he the ordering 
of Q in the counterclockwise direction, see Fig. 3. Noting that the real part Rc{vmi) 
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can be written in terms of the dot product, Re{vmi) = v^'h (2.24) can be recast as 

J2 9™ e''^-(^)+<-3 = 0, (2.30) 

where (fy'^iL) = <f)m{L) + lm{vmi). 

On the basis of (2.30), it is easy to verify the following facts: Let 3 = j^je, with e a 
unit vector in C. An inspection of (2.30) shows that, for I3 1 ^ 1, the roots of (2.30) will 
concentrate along the "directions" for which the projection of e on at least two w* 's is 
the same. Invoking the convexity assumption (Assumption A4), this can only happen 
when ?;* • e = • e for some n. In such cases, the contributions of the terms with 
indices m 7^ n, n + 1 in (2.30) are negligible — at least once \^\ ^ 1 — and the zeros 
will thus asymptotically lie along the half-lines given in the parametric form by 

3 = 3ft) = ij^rj"^ 12 M^) + ^ti< - *e[0,oo). (2.31) 

Clearly, the latter is a line perpendicular to the (n, n + l)-st side of the convex set with 
vertices vl, . . . ,v*, which is shifted (away from the origin) along the corresponding 
side by a factor proportional to log(g„+i/(7„), see Fig. 3. 

Sufficiently far away from zm (on the scale L~'^), the zeros resume the pattern es- 
tablished around the two-phase coexistence curves. In particular, the zeros are asymp- 
totically equally spaced but their overall shift along the asymptote is determined by the 
factor (prn{L) — which we note depends very sensitively on L. Computer simulations 
show that, at least in generic cases, this pattern will persists all the way down to the 
multiple point. Thus, even on the "microscopic" level, the zeros seem to form a "phase 
diagram." However, due to the lateral shifts caused by qm+i 7^ 3m. a "macroscopic" 
quadruple point may resolve into two "microscopic" triple points, and similarly for 
higher-order multiple points. 



3. Characterization of phase diagrams 

The goal of this section is to give the proof of Theorem 2.1. We begin by proving a series 
of auxiliary lemmas whose purpose is to elevate the pointwise Assumptions A3-A4 into 
statements extending over a small neighborhood of each coexistence point. 



3.1. Auxiliary claims. Recall the definitions of S^m, Q{z) and Vm{z), in (2.2), (2.5) 
and (2.23), respectively. The first lemma gives a limiting characterization of stabiUty of 
phases around coexistence points. 

Lemma 3.1. Let Assumption A1-A2 hold and let z G & be such that |Q(z)| > 2. 
Let (z-k) he a sequence of numbers Zk & & such that Zk ^ z but Zk ^ z for all k. 
Suppose that 

= lim (3.1) 

fe— >c» \Zk — Z\ 

exists and let m G Q{z). If Zk € S^mfor infinitely many fc > 1, then 

Re(e'%„) > Re(e'^z;„) for all n G Q{z), (3.2) 
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wherevn = Vn{z). Conversely, if the inequality in(3. 2) fails for at least one n G Q{z), 
then there is ane > such that 

W,fi{z) = {zGe:\z-z\<^,zi^z, l^-e'^l < e} (3.3) 

has empty intersection with J^m, i.e., H W^fii^^ = 0. In particular, Zk ^ S^mfor k 
large enough. 

Remark 5. In the following, it will be useful to recall some simple facts about complex 
functions. Let /, g and h be functions C ^ C and let dz and dg be as in (2.1). If / 
satisfies dzf{zo) = (i.e., Cauchy-Riemann conditions), then all directional derivatives 
of / at 2;o = xo+ iyo can be expressed using one complex number A = dzf{xo + iyo). 
i.e., we have 

/ (.To + ecos(p + ij/o + iesinc/j) - /(.To + iyo) = e-4e"^ + o(e), e 4 0, (3.4) 

holds for every (p G [— tt, tt). Moreover, if g is differentiable with respect to x and y at 
zq = xq + iyo and h satisfies dzh{z') = at z' = g{zo), then the chain rule holds for 
z h{g{z)) atz = zq. In particular, dzh{g{zo)) = {dzh){g{zo))dzg{zo). 

Proof of Lemma 3.1. Let m e Q{z) be fixed. Whenever Zk € J^m. we have 

l0g|Cm(2;fe)| - l0g|Cm(^)| > log|Cn(^fe)| " log|Cn(^)|, « € Q{z), (3.5) 

because |Cm(^)| = |Cn(^)l> by our assumption that 171,71 G Q{z). Using the notation 

= ^ (3.6, 

for n e Q{z) (which is well defined and non-zero in a neighborhood of z), the inequal- 
ity (3.5) becomes 

log|F„,„(zfc)| - log|F„,„(^)| > 0, n e Qiz). (3.7) 

Note that the complex derivative dzFm,n{z) exists for all n e Q{z). Our task is then to 
prove that 

Re(e'« y-'-^f ) > 0, n e Q{z). (3.8) 

Fix n G Qiz). Viewing z ^ Fm.n{z) as a function of two real variables x = Rez and 
y — Imz, we can expand log \Fm,n{z) \ into a Taylor series around the point z to get 

\og\FmM\ - log|F„,„(0)| = Ref (zfe - z)^^^^^^) + 0{\zk - zf). (3.9) 

V ^m,n[Z) J 

To derive (3.9) we recalled that Fm^n is at least twice continuously differentiable (hence 
the error bound) and then applied the identity 

— d-x — + — d-y — ""y^ = ""v^^ - '^^^ 

where Axk = Re(zfe — z) and Ayk = Im(zfe — z). (To derive (3.10), we just have to 
apply the chain rule to the functions z 1-^ log Fm,n{z). See Remark 5 for a discussion 



(3.10) 
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of this point.) Using that Zk z, the inequality (3.8) and hence also (3.2) now follows 
by combining (3.9) with (3.5), dividing by \zk — z\ and taking the limit k ^ oo. 

If, on the contrary, the inequahty (3.2) is violated for some n e then (3.8) 

fails to hold as well and hence (3.7) and (3.5), with Zk replaced by z, must be wrong for 
z e We^0{z) whenever e is small enough. But m e Q{z) implies that |Cm(^)| = |Cn(^)| 
and thus |Cm(2;)| < \Cn{z)\ for all 2; e W^fii^z), proving that ^™ n W^fi{z) = 0. By 
(3.1) and the fact that Zk — » z, we have Zk € W^fi{z) and hence Zk for all k 
large enough. □ 

Lemma 3.1 directly implies the following corollary. 

Corollary 3.2. Let Assumption A1-A2 hold and let m,n € TZ be distinct. Let {zk) be 
a sequence of numbers Zk € H S^n such that Zk ^ z & but Zk 7^ zfor all k. 
Suppose that the limit (3.1) exists and equals e'^. Then ^^{e^^Vm) = Re(e'^w„). 

Proof. Follows iimnediately applying (3.2) twice. □ 

The next lenmia will ensure that multiple points do not cluster and that the coexis- 
tence lines always intersect at positive angles. 

Lemma 3.3. Suppose that Assumption A holds and let z &. Suppose there are two 
sequences {zk) and {z'l^) of numbers from such that \zk — ^1 = — z| 7^ for 
all k and Zk, z'f. ^ z as k ^ oo. Let a,b, c G TZ and suppose that zt e J>^a H and 
z'j. e n S^cfor all k. Suppose the limit (3.1) exists for both sequences and let e'^ 
and e'^ be the corresponding limiting values. 

(1) Ifa,b,c are distinct, then e'^ 7^ e'^ . 

(2) If a ^ h = c and Zk 9^ z'j. for infinitely many k, then \ Q{z) \ = 2 and e'^ = — e'^ . 

Remark 6. The conclusions of part (2) have a very natural interpretation. Indeed, in 
this case, z is a point on a two-phase coexistence line (whose existence we have not 
estabUshed yet) and Zk and z'j. are the (eventually unique) intersections of this Une with 
a circle of radius \zk — z\ = \z'i, — z\ around z. As the radius of this circle decreases, 
the intersections Zk and z^ approach z from "opposite" sides, which explains why we 
should expect to have e'^ = — e'^ . 

Proof of Lemma 3.3. Throughout the proof, we set Vm = Vm{z). We begin by prov- 
ing (1). Assume that a,b,c e 7?. are distinct and suppose that e'^ = e'^ . Note that, 
since Q{z) D {a, h, c}, the point ^ is a multiple point. Corollary 3.2 then implies that 

Re(e'%a) = Re(e'%6) = Re(e'%c), (3.11) 

and hence Va, Vb and Vc lie on a straight line in C. But then Vb and Vc cannot simulta- 
neously be vertices of a strictly convex polygon, in contradiction with Assumption A4. 

In order to prove part (2), let a 7^ = c, suppose without loss of generality that 
Zk 4 for all k. If e'^ ^ ±e'^', then Corollary 3.2 implies that Re(e'''(-ya - Vb)) = 
= Re(e'^ {va — Vb)) and hence Va = Vb, in contradiction with Assumption A3. Next 
we will rule out the possibility that e'^ = e'^ , regardless of how many phases are stable 
at z. Let G{z) = Ca{z)/Cb{z) and note that \G{zk)\ = 1 = |G(4)I for all k. Applying 
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Taylor's theorem (analogously to the derivation of (3.9)), dividing by — and 
passing to the hmit fc — > oo, we derive 

lin,Ref^^%g(£^Uo. (3.12) 



k^oo \\zk-z'f,\ G{zk) 

The second ratio on the left-hand side tends to Va — Vb. As for the first ratio, an easy 
computation reveals that, since \zk — z\ = \z'j^ — z\ ^ Q, we have 

^H{0k-e'^)/2) ^ ^^^^^ 



|sin((efe-W)l 
where 

g.e. = and e< = ^f^. (3.14) 

\zk - z\ \z'^ - z\ 

By our assumptions, we have e'®'= e'^ and e'^'' ^ e'^ as fc ^ oo. Suppose now that 
gi0 _ gi0 xhen, choosing a subsequence if necessary, the left-hand side of (3.13) tends 
to a definite sign times ie'^. Inserting this into (3.12) and using Corollary 3.2, in addition 
toRe(e'^(ua — Wfe)) = 0, wenowgetthatalsoRe(ie'^(i;o — 1^5)) = Im(e'^ (ua — Ub)) = 0. 
Consequently, Va = Vb, again contradicting Assumption A3. 

To finish the proof of the claim (2), it remains to rule out the possibility that e'^ = 
— e'^ in the case when is a multiple point. Let n G Q{z) be another phase stable at z, 
i.e., a,b. By Lemma 3.1, we have 

Re(e'^(u„ - w„)) > and Re(e'^'(i;„ - w„)) > 0, m = a,b. (3.15) 

But then e'^ = — e'^ would imply that Re{e^^Va) = Re(e'^i;„) = Re(e'^i;b), in contra- 
diction with Assumption A4. Therefore, | Q(^) | < 3, as claimed. □ 

Corollary 3.4. Suppose that Assumption A holds and let z € he a multiple point. 
Then there exists a constant S > such that \Q{z)\ < 2 for all z € {z' G 0': < 
\z' — z\ < 6}. In particular, each multiple point in G is isolated. 

Proof. Suppose z G ^ is a non-isolated multiple point. Then there is a sequence z^G G 
such that Zh ^ z and, without loss of generality, Q(zfe) = Qo with \Q.(i\>2>, z^ z 
for all k, and such that the limit (3.1) exists. Taking for (z^) the identical sequence, 
z'l^ = Zk, we get e'^ = e'^ in contradiction to Lermna 3.3(1). Therefore, every multiple 
point in is isolated. □ 

Our last auxiliary claim concerns the connectivity of sets of 6 such that (3.2) holds. 
As will be seen in the proof of Lemma 3.6, this will be crucial for characterizing the 
topology of the phase diagram in small neighborhoods of multiple points. 

Lemma 3.5. Suppose that Assumption A holds and let z £ G be a multiple point. For 
m e Q{z), let Vm = Vm{z). Then, for each m G Q{z), the set 

Im = {e'^ : e e [0, 2tt), Re{e'%m) > Re(e'^t;„), n e Q(^) \ {m}} (3.16) 

is connectedand open as a subset of{z G 0: \z\ = 1}. In particular, ife^^ is such that 

Re{e'^Vm) ^ max Re(e'^t;„), (3.17) 

then e'^ is one of the two boundary points oflm- 
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Proof. By Assumption A4, the numbers m G Q{z), are the vertices of a strictly 
convex polygon ^ in C. Let s = \Q{z)\ and let (fi, . . . , Vs) be an ordering of the 
vertices of ^ in the counterclockwise direction. For m = 1 , . . . , s define Avm = Vm — 
Vm-i, where we take va = Note that, by strict convexity of the arguments 9„i 
of Avm, i.e., numbers 0m such that Avm = \Avm\e^^'", are such that the vectors 
e'^i, . . . , e^^" are ordered counterclockwise, with the angle between e'^"" and e'^^+i 
lying strictly between and tt for all m = 1, . . . ,s (again, we identify m = 1 and 
m = s + 1). In other words, for each to, the angles 6i ... ,6s can be chosen in such a 
way that 9m < Om+i < ■ ■■ < 6m+s, with < 9m+k - Om+k-l <'K,k = l,...,S. 
(Again, we identified m + k with m + k — s whenever m + k > .s). 

Using Jm to denote the set Jm = {ie~"' : i? S {9m, dm+i)}, we claim that Im = 
Jm for all TO = 1, . . . , s. First, let us show that Jm C Im- Let thus i? e {9m, Om+i) and 
observe that 

Re{ie-'^ Avm) = \Avm\sm{'& - 9m) > 0, (3.18) 
because 9m < Om+i < 9m + '^- Similarly, 

Re(ie-'''Zit;„+i) = \Avm+i\ sin{i) - 9m+i) < 0, (3.19) 

because 9m+i — tt < 0m < ^ < 0m+i- Consequently, Re{ie~^'^Vm) > Re(ie~"'i'„) 
holds for both n = m + 1 and n = m — 1. 

It remains to show that Re{ie~^'^Vm) > Re(ie~"'ii„) is true also for all remaining 
n G Q{z). Let n G Q{z) \ {m, m ± 1}. We will separately analyze the cases with 
0n — 0m G (0,7r] and 6'„ — e (— 7r,0). Suppose first that 6'„ — 0m € (0, tt]. This 
allows us to write n = to + for some k e {2,...,s — 1} and estimate 

k 

Re(ie-'''K - Vm)) = ^Re(ie-'''zit;„+,) 

k 

= \AVm+MH^ - 0m+j) < 0. (3.20) 

The inequality holds since, in light of < 0m+i < ■ ■ ■ < 0m+k < + tt, each sine is 
negative except perhaps for the last one which is allowed to be zero. On the other hand, 
if 0n — 0m & {~'^, 0)' we write n = m — k instead, for some k & {2, . . . , s — 1}, and 
estimate 



Re(ie-"'K-t;„))= ^ Re(ie-"'/A^;„+,■) 
j=-k+i 



= \^Vm+j\sm{^-0m+j)>O- (3.21) 

j=-k+l 

Here we invoked the inequalities — tt < 0m -k < ■ ■ ■ < 0m < to show that each 
sine on the right-hand side is strictly positive. 

As a consequence of the previous estimates, we conclude that Jm C Im for all 
TO = 1, . . . , s. However, the union of all Jm's covers the unit circle with the exception 

of s points and, since the sets Im are open and disjoint, we must have /,„ = Jm for 
all TO G Q{z). Then, necessarily, Im is connected and open. Now the left-hand side of 
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(3.17) is strictly greater than the right-hand side for e G /,„, and strictly smaller than 
the right-hand side for e'^ in the interior of the complement of /„• By continuity of 
both sides, (3.17) can hold only on the boundary of . □ 



3.2. Proof of Theorem 2.1. Having all the necessary tools ready, we can start prov- 
ing Theorem 2.1. First we will apply Lemma 3.5 to characterize the situation around 

multiple points. 

Lemma 3.6. Suppose that Assumption A holds and let z & G he a multiple point. For 
5 > 0, let 

= 1^ g ^. 1^ _ ^1 = Q(^) 3 (3_22) 
Then the following is true once 5 is sufficiently small: 

(1) For each m G Q{z), the set is connected and has a non-empty interior 

(2) I'ir} = whenever m ^ Q{z). 

(3) For distinct m and n, the sets Im^ and In^ intersect in at most one point. 

Proof. The fact that Im'' = for m ^ Q{z) once (5 > is sufficiently small is a 
direct consequence of the continuity of the functions (^m and Indeed, if there were a 
sequence of points Zk tending to z such that a phase m were stable at each Zk, then m 
would be also stable at z. 

We will proceed by proving that, as ^ J. 0, each set Im\ m G Qiz), will even- 
tually have a non-empty interior. Let m G Q{z). Observe that, by Lemma 3.5, there 
is a value e'^ (namely, a number from /„) such that Re(e'^Wm) > Re(e'^w„) for all 
n G Q{z) \ {m}. But then the second part of Lemma 3.1 guarantees the existence of 
an e > such that Q{z) = {m} for all z G We^0{z) — see (3.3). In particular, the in- 
tersection We^0{z) r\ {z: \z — z\ = 6}, which is non-empty and (relatively) open for 

(5 < e, is a subset of Im^ . It follows that the set Im^ has a nonempty interior once S is 
sufficiently small. 

Next we will prove that each Im, m G Q(z), is eventually connected. Suppose 

that there exist a phase a G Qiz) and a sequence i5fc i such that all sets 1^^'°^ are 

not connected. Then, using the fact that la''^ has nonempty interior and thus cannot 
consist of just two separated points, we conclude that the phase a coexists with some 
other phase at at least three distinct points on each circle {z: \z — z\ = 6k}. Exphcitly, 

there exist (not necessarily distinct) indices b^j^^ G Q{z) \ {a} and points (z'^''), j = 
1,2,3, with \z[^^ - z\ = 6k and 4^'^ -^f^ for j ^ £, such that a,b\^^ G 2(4^^)- 
Moreover, (choosing subsequences if needed) we can assume that b[^^ for some 

6^^-' G Q{z) \ {a} independent of k. Resorting again to subsequences, we also may 
assume that the limits in (3.1) exist for all three sequences. 

Let us use e'^^ to denote the corresponding hmits for the three sequences. First we 
claim that the numbers e'^^ , j = 1, 2, 3, are necessarily all distinct. Indeed, suppose 
two of the e'^' 's are the same and let b and c be the phases coexisting with a along 
the corresponding sequences. Then Lemma 3.3(1) forces b = c, which contradicts both 
conclusions of Lemma 3.3(2). Therefore, all three e'^^ must be different. Applying now 
Corollary 3.2 and Lemma 3.1, we get Re(e'®JWa) = max„gQ(2)\{(j} Re(e'^^^„) for 
j = 1, 2, 3. According to Lemma 3.5, all three distinct numbers e'^^ , j = 1, 2, 3, are 
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endpoints of /„, which is not possible since /„ is a connected subset of the unit circle. 

Thus, we can conclude that /i*^ must be connected once (5 > is sufficiently small. 

To finish the proof, we need to show that la ^ fl I^^^ contains at most one point for 
any a ^ b. First note that we just ruled out the possibility that this intersection contains 
three distinct points for a sequence of S's tending to zero. (Indeed, then a would coexist 
with b along three distinct sequences, which would in turn imply that a and b coexists 
along three distinct directions, in contradiction with Lemma 3.5.) Suppose now that 
la^^ n ij^^^ contains two distinct points. Since both li^^ and 4'') are connected with 
open interior, this would mean that ^ and cover the entire circle of radius S. 
Once again, applying the fact that two Im have at most two points in common, we 
then must have li^^ = for all c ^ a, 6. But Q{z) contains at least three phases which 

necessitates that ■* ^ for at least three distinct m. Hence la^ fl ijf^ cannot contain 
more than one point. □ 

Next we will give a local characterization of two-phase coexistence lines. 

Lemma 3.7. Suppose that Assumption A holds and let m,n TZ be distinct. Let z 
he such that z G S^m H S^n and Q,(z') C {TO,n} for z' G ©^(z). Then there exist 
numbers 5' £ (0,(5), ti <Q,t2> 0, and an twice continuously dijferentiable function 
7z : (^1,^2) ^5'{z) such that 

(1) 7z(0) = z. 

(2) Kmii^m = Knh.m = Chzit)), t e {hM 

(3) \imtit, lz{t), limt|t2 72 (t) G dns'{z). 

The curve t t-^ (t) is unique up to reparametrization. Moreover, the set D^' {z) \ 
lz{ti 1 ^2) has two connected components and m is the only stable phase in one of the 
components while n is the only stable phase in the other. 

Proof. We begin by observing that by Assumption A3, the function 

4>m,n{x,y) = log|Cm(a; + i?/)| -log|Cn(a; + iy)| =Ke,\ogFm,n{x + iy)., (3.23) 

has at least one of the derivatives dxd>m.mdy4>m,n non-vanishing at x + iy = z. 
By continuity, there exists a constant 77 > such that one of the derivatives is uni- 
formly bounded away from zero for all z' = u + w e ^•qiz). Since z = x + \y £ 
,5^„i n ,5^n, we have (f>m.n{x,y) = 0. By the implicit function theorem, there exist 
numbers t[^, t[, xq, xi, yo and t/i such that < < t'l, xq < x < xi, yo < y < 
yi and (.xo,.ti) x (2/0,2/1) C IQ',,(z), and twice continuously differentiable functions 
u: (to.i'i) (xq,Xi) andw: (i'o.fi) {yo,yi) such that 

</'m,nKt),t^(i)) =0, te{t'o,t[), (3.24) 

and 

u{0) ^ X, and v(0) = y. (3.25) 

Moreover, since the second derivatives of 0m, n are continuous in if? and therefore 
bounded in D,, (z), standard theorems on uniqueness of the solutions of ODEs guarantee 
that the solution to (3.24) and (3.25) is unique up to reparametrization. The construction 
of jz is now finished by picking 6' so small that D^/ (z) c (xq , xi ) x (yo , 2/1 )> and tak- 
ing to and ti to be the first backward and forward time, respectively, when {u{t),v{t)) 
leaves B)g/{z). 
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The fact that IS>s' (z) \ ■jz {ti , ^2) splits into two components is a consequence of the 
construction of 7^. Moreover, 7^ is a (zero-)level curve of function (j)m,n which has a 
non-zero gradient. Hence, 0m, „ < on one component of Hg'iz) \ 7z(ii, t2), while 
0m, n > on the other. Recalling the assumption that Q{z') C {m,n} for z' in a 
neighborhood of z, the claim follows. □ 

Now we can finally give the proof of Theorem 2.1. 
Proof of Theorem 2.1. Let ^ denote the set of all multiple points in ^, i.e., let 

= {zeff: \Qiz)\ > 3}. (3.26) 

By Corollary 3.4, we know that ^ is relatively closed in and so the set 0' = &\ ^ 
is open. Moreover, the set 'S ^ G' consists solely of points where exactly two phases 
coexist. Lemma 3.7 then shows that for each 2; G fl there exists a disc D^/ {£) and 
a unique, smooth in '^(,i{z) passing through z such that Q[z') = Q,{z) for all z' on 
the curve 7^. Let 7^ be a maximal extension of the curve 7,, in 0' . We claim that 7^ is 
either a closed curve or an arc with both endpoints on 9^'. Indeed, if 7^ were open with 
an end-point z G 6' , then Q.[z) D Q,{z), by continuity of functions Cm- But z G G' 
and so |Q(z)| < 2, which implies that Q,{z) = Q{z). By Lemma 3.7, there exists a 
non-trivial curve 75 along which the two phases from Q{z) coexist in a neighborhood 
of z. But then 7^ U would be a non-trivial extension of -jz, in contradiction with the 
maximality of 72. Thus we can conclude that z G dG". 

Let e denote the set of maximal extensions of the curves {72 : 2 e $f fl ^' } . Let ^ C 
^ be a compact set and note that Corollary 3.4 implies that is finite. Let Sq be so 

small that, for each 2m € ^fl^, wehaveD5(,(zM) C ^,D5o(zM)ri^ = {^m} and the 
statements in Lemma 3.6 hold true for 6 < So- Let S € (0, i5o]. We claim that if a curve 
"^/f E Q intersects the disc ©^(zm) for ^ -^m G H i^, then the restriction '^/f n ©^(zm) 
is a simple curve connecting zm to ^©^(zm)- Indeed, each curve 'rf G G terminates 
either on dG or on If "enters" IDsizu) and does not hit zm, our assumptions 
about 6a imply that "leaves" B>s{zm) through the boundary. But Lemma 3.7 ensures 
that one of the phases coexisting along dominates in a small neighborhood on the 
"left" of while the other dominates in a small neighborhood on the "right" of 
The only way this can be made consistent with the connectivity of the sets /m^ in 
Lemma 3.6 is by assuming that Im^ ^ only for the two m's coexisting along 'rf. But 
that still contradicts Lemma 3.6, by which 7m' ^ for at least three distinct m. Thus, 
once a curve G Q intersects ©^(zm), it must terminate at zm- 

Let ^0 = \ Uze^ ^So (z) and let Z\ : ^0 ^ [0, 00) be a function given by 

A{z) = M{6' e (0, 5o) : ©5/ {z) C G, ©5/ (z) n U^f ee is disconnected}. (3.27) 

We claim that A is bounded from below by a positive constant. Indeed, A is clearly 
continuous and, since &o is compact, A attains its minimum at some z £ !^q. If 
A{z) = 0, then 2; is a limit point of [J<^^q and thus z G^ for some € C. More- 
over, for infinitely many 6' G (0, (5o), the circle d^f,i(z) intersects the set U-rTee 
in at least three different points. Indeed, the curve "if 9 z provides two intersections; 
the third intersection is obtained by adjusting the radius 6' so that ©5'(z) n U'^f ee 
is disconnected. Thus, we are (again) able to construct three sequences (z^), (zj^) 
and (z^') such that, without loss of generahty, z^, z^, z'l G 5^a H 5^h for some dis- 
tinct a,b G TZ (only two phases can exist in sufficiently small neighborhoods of points 
in \zk-z\ = |z^ -z\ = \z'l - z| 0, but z^^ z'^^ z'{. ^ z^ for all k. However, 
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this contradicts Lemma 3.3, because its part (2) cannot hold simultaneously for all three 
pairs of sequences {zk, z'f.), {z'l^, z'l) and (zfc, z'l). 

Now we are ready to define the set of points zi , . . . , z^. Let e be the minimum of the 
function A in SIq and let 5 = niin((5o, e). Consider the following collections of open 
finite discs: 



It is easy to check that the union of these discs covers (i). Let 8 = §i U §2 U S3. By 
compactness of '3, we can choose a finite collection §' c § still covering 3. It remains 
to show that the sets ^ = n D for D e §' will have the desired properties. Let D e S' 
and let z be the center of ©. If ID € S3, then n © = 0. Indeed, if z' is a coexistence 
point, then ©^(z') e Si U S2 and thus dist(z, z') > (5 + \b and hence z' ^ O. Next, 
if D G S2, then z G J# and, by the definition of 5q and e, the disc D contains no multiple 
point and intersects ^ only in one component. This component is necessarily part of 
one of the curves e C. Finally, if ED e §1, then z is a multiple point and, relying on 
our previous reasoning, several curves € C connect z to the boundary of U. Since 
Lemma 3.6 imphes the existence of exactly | Q{z) \ coexistence points on 9D, there are 
exactly such curves. The proof is finished by noting that every multiple point 

appears as the center of some disc P e S', because that is how the collections (3.28) 
were constructed. □ 



4. Partition function zeros 

The goal of this section is to prove Theorems 2.2-2.5. The principal tool which en- 
ables us to control the distance between the roots of Z^^^ and the solutions of equations 
(2.17-2.18) or (2.24) is Rouche's Theorem (see e.g. [16]). For reader's convenience, we 
transcribe the corresponding statement here: 

Theorem 4.1 (Rouche's Theorem). Let !3 C C be a bounded domain with piecewise 
smooth boundary 82!. Let f and g be analytic on2UdS>.If\g{z)\ < \f{z)\ for all z € 
d2>, then f and f + g have the same number of zeros in 2, counting multiplicities. 

More details on the use of this theorem and the corresponding bounds are stated in 
Sect. 4.2 for the case of two-phase coexistence and in Sect. 4.4 for the case of multiple 
phase coexistence. 

Root degeneracy will be controlled using a link between the non-degeneracy condi- 
tions from Assumption B and certain Vandermonde determinants; cf Sect. 4.1. Through- 
out this section, we will use the shorthand 



§1 
S2 
§3 



{^f,{z): z e 

{©^(z) : z e n U^ge dist(z, Ud^s, ^) > ¥} 
{Bs{z):ze2, dist(z,UDe8iUS. ») > 1^}- 



(3.28) 




(4.1) 



meQ 



to denote the set of points z G ^ where all phases from a non-empty Q C TZ are 
"almost stable" (as quantified by e > 0). 
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4.1. Root degeneracy. In this section we will prove Theorem 2.2. We begin with a 
claim about the Vandermonde matrix defined in terms of the functions 

hm{z) = ^4i^^ ^^'^^/dml (4.2) 

where the dependence of bm on L has been suppressed in the notation. Let us fix a 
non-empty Q C TZ and let q = \Q\. For each z e ^k/l(Q), we introduce the g x g 
Vandermonde matrix M(^;) with elements 

Me,m{z) = bm{zY, me Q, £ = 0,l,...,g-l. (4.3) 

Let ||M|| denote the £^(Q)-norm of M (again without making the Q-dependence of this 
norm notationally expUcit). ExpUcitly, ||M|p is defined by the supremum 



: V |w„|2 = 1 I (4.4) 



q-l 

e=o meQ meQ 

where (wm) is a |Q| -dimensional complex vector. 

Throughout the rest of this section, the symbol || • j| will refer to the (vector or 
matrix) i^-^-norm as specified above. The only exceptions are the £P-norms ||q||i, ||q||2 
and ||q||oo of the r-tuple ((j'm)meK. which are defined in the usual way. 

Lemma 4.2. Suppose that Assumption B3 holds and let Lq be as in Assumption B3. For 
each Q CTZ, there exists a constant K = K{Q) < oo such that 

\\M-\z)\\ < K, for all ze y^/L{Q) andL> Lq. (4.5) 

In particular, yi{z) is invertible for all z G =5^k/l(Q) cind L > Lq. 

Proof Let Q C 7?. and g = Let us choose a point z G =5^k/l(2) and let M 
and hm, me Q, be the quantities M(2;) and hm{z), m G Q. First we note that, 
since M is a Vandermonde matrix, its determinant can be explicitly computed: dct M = 
UnKni^n — ^m)> where "<" dcuotcs a complete order on Q. In particular. Assump- 
tion B3 impUes that |detM| > a«(«-i)/2 > once L > Lq. 

To estimate the matrix norm of M~^, let Ai , . . . , be the eigenvalues of the Her- 
mitian matrix MM+ and note that A^ > for all ^ — 1, . . . , g by our lower bound 
on |dctM|. Now, ||M+jp is equal to the spectral radius of the operator MM+, and 
||M~^|p is equal to the spectral radius of the operator (MM+)~^. By the well-known 
properties of the norm we thus have 

||Mf = ||M+f = max A^, (4.6) 



- max XJ^. (4.7) 



while 

||M-i||2 = 

Now |det = det M M+ = Ai . . . Ag and a simple algebraic argument gives us that 
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Using the lower bound on |detM|, this implies that \\M-^\\ < V- ||M||9"^ The 
claim then follows by invoking the uniform boundedness of the matrix elements of M 
(see the upper bound from Assumption B3), which impUes that ||M|| and hence also 
||M~^ II is uniformly bounded from above throughout ^k/l{Q)- O 

Now we are ready to prove Theorem 2.2. To make the reading easier, let us note that 
for Q = {to}, the expression (2.8) defining '^e(Q) can be simplified to 

%{{m}) = {z€^: |C„(^)| <e-^/2|C(^)|foraUn^m}, (4.9) 
a fact already mentioned right after (2.8). 

Proof of Theorem 2.2. Let m G TZ. Since the sets '^k/l (Q)^ Q C TZ, cover i?, it suffices 
to prove that 7^ in '^l-^u^lH'^}) '?^n/L{Q) for each Q C 7^. In fact, since 
z € ^L-''a)t ({"^l) implies that to is stable, |(^„i(z)| = C(^)' we may assume without 
loss of generality that m G Q, because otherwise '^i-<i(^^({TO})n'^K/j;,(Q) = 0.Thus, 
let to e Q c 7?. and fix a point z € '^l-'ul ({"^}) ^ '^k,/l{Q)- By Assumption B4, 
we have the bound 



\zr{z)\>az)' * 



- ^ Qn 

n£C\{m} 



-CoL'^llqllie- 



tL 



(4.10) 



Since z G ^l—^ul ({"^})' ^^^^ ICn(-2^)| < C(^)e ''"^ for n ^ to. In conjunc- 
tion with Assumption B2, this impUes 



Ciz) 



tH-tL 1, 



n^m. 



(4.11) 



On the other hand, we also have 



6^\zi 



> e 



(4.12) 



where we used that !Crn(2:)| = (^{z). Since cjl x, (4.11^.12) show that the right- 
hand side (4.10) is dominated by the term with index to, which is bounded away from 
zero uniformly in L. Consequently, Z^'^ 7^ throughout '^i,-<ia;i,({™}) ^ '^k/l(Q). 
provided L is sufficiently large. 

Next we will prove the claim about the degeneracy of the roots. Let us fix Q c 7?. 
and let, as before, q = \ Q\. Suppose that L> Lq and let 2 € ^^^/l ( Q) be a root of Z^'^ 
that is at least g'-times degenerate. Since Z^'^ is analytic in a neighborhood of z, we 
have 

diZf\z) = 0, i = 0,l,...,q-l. (4.13) 

It will be convenient to introduce g-dimensional vectors x = x{z) and y = y{z) such 
that (4.13) can be expressed as 

M(2)x = y, (4.14) 
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with M(^;) given by (4.2) and (4.3). Indeed, let x = x(^) be the vector with components 

Xm = qJ^^^rV-) , meQ. (4.15) 
Similarly, let y = y(z) be the vector with components yo, ■ ■ ■ , Uq-i, where 

(4.16) 

Recalling the definition ^q^l{z) from (2.13), it is easily seen that (4.14) is equivalent 
to (4.13). 

We will now produce appropriate bounds on the -norms ||y|| and ||x|| which 
hold uniformly in z G "^k/lCQ)^ and show that (4.14) contradicts Lemma 4.2. To 
estimate 1 1 y j | , we first note that there is a constant ^ < oo, independent of L, such that, 
for alH = 0, . . . , g - 1 and aU ^ e '^«/l(Q), 

L-^%[c^j:\z)f - hm{zY[c^J:\z)f < AL-'az)^\ {aid 

Here the leading order term from L~'^^'dl[Cm\z)Y' is exactly canceled by the term 
bm{zY[C,m\z)]^'^ , and the remaining terms can be bounded using (2.11). Invoking 
(4.17) in (4.16) and applying (2.14), we get 

||y|| < AM^^L-^' + max_^Q)||q||lV«iV^^ (4.18) 

where the factor ^ comes from the conversion of ^^-type bounds (4.17) into a bound 
on the £^-norm ||y||. On the other hand, by (2.9) and gm > 1 we immediately have 

||x|| > e-^"""". (4.19) 

But ||x|| < ||M^^(z)|| ||yj|, so once L is sufficiently large, this contradicts the upper 
bound ||M^^(z)|| < if impUed by Lemma 4.2. Therefore, the root at ^; cannot be more 
than (g — l)-times degenerate after all. □ 



4.2. Two-phase coexistence. Here we will prove Theorem 2.3 on the location of parti- 
tion function zeros in the range of parameter z where only two phases from 1Z prevail. 
Throughout this section we will assume that Assumptions A and B are satisfied and use 
K and r to denote the constants from Assumption B. We will also use 5l{z) for the 
function defined in (2.16). 

The proof of Theorem 2.3 is based directly on three technical lemmas, namely. 
Lemma 4.3-4.5 below, whose proofs are deferred to Sect. 5.2. The general strategy 
is as follows: First, by Lemma 4.3, we will know that the solutions to (2.17-2.18) 
are within an 0(e^'"^)-neighborhood from the solutions of similar equations, where 

the functions Cm get replaced by their analytic counterparts C,m \ Focusing on specific 
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indices m and n, we will write these analytic versions of (2.17-2.18) as f{z) = 0, 
where / is the function defined by 

f{z)=qmC^\z)''' +qnCL''H^)''\ ze.y^/Li{m,n}). (4.20) 

The crux of the proof of Theorem 2.3 is then to show that the solutions of f{z) = are 
located within an appropriate distance from the zeros of Z^"{z). This will be achieved 
by invoking Rouche's Theorem for the functions / and f + g, where g is defined by 

g{z) = Zf\z) - f{z), z e y^/L{{m, n}). (4.21) 

To apply Rouche's Theorem, we will need that \g{z) \ < \f{z) \ on boundaries of certain 
discs in n}); this assumption will be verified by combining Lemma 4.4 (a 

lower bound on |/(z)|) with Lemma 4.5 (an upper bound on [(/(-z)]). The argument 
is then finished by applying Lemma 4.3 once again to conclude that any two distinct 
solutions of the equations (2.17-2.18), and thus also any two distinct roots of Z^^, are 
farther than a uniformly-positive constant times L"''. The actual proof follows a sUghtly 
different path than indicated here in order to address certain technical details. 

We begin by stating the aforementioned technical lemmas. The first lemma provides 
the necessary control over the distance between the solutions of (2.17-2.18) and those 
of the equation f{z) = 0. The function / is analytic and it thus makes sense to consider 
the multiplicity of the solutions. For that reason we will prefer to talk about the roots of 
the function /. 

Lemma 4.3. There exist finite, positive constants Bi, B2, Ci and Li, satisfying the 
bounds El < B2 and Cie^"^^ < BiL^'^ whenever L > Li, such that for all L > L\, 
all s < {Bi + B2)L^'^ and all G yK/{2L){{m,n}) with IDs(-Zo) C the disc 
Ds(zo) is a subset of ^^/[^{{m, n}) and the following statements hold: 

(1) If s < BiL~''^, then disc 10^(2:0) contains at most one solution of the equations 
(2.17-2.18) andatmost one root of function f, which is therefore non-degenerate. 

(2) If s > Cie"^^ and if zq is a solution of the equations (2.17-2.18), then ©^(2:0) 
contains at least one root of f. 

(3) Ifs > C\e~'^^ and if zq is a root of the function f, then ^s{zo) contains at least 
one solution of the equations (2.17-2.18). 

(4) Ifs = B2L~'^ and if both m and n are stable at zq, then ©^(^o) contains at least 
one solution of the equations (2.17-2.18). 

The next two lemmas state bounds on |/(z)| and \g{z)\ that will be needed to apply 
Rouche's Theorem. First we state a lower bound on |/(z) |: 

Lemma 4.4. There exist finite, positive constants 62 and C2 obeying c.2 < C2 and, for 
any C > C2 and any sequence (e^) of positive numbers satisfying 

lim L'^CL = 0, (4.22) 

L — >oo 

there exists a constant L2 < 00 such that for all L > L2 the following is true: If zq 
is a point in S^kK^l) ({^n, n}) fl [S^m U <5^n) and ©c^^ (^o) C G, then there exists a 
number s{zo) e {c2eL,C2eL} suchthatI]>g(^;^g){zo) C S^^/ {2 l) {{1^,71}) and 

liminf inf > e^i'^CC^o)^'- (4.23) 

sTs(zo) z: \z-zo\=s 

Moreover, if f has a root in ©026^ (-^o)- then s{zo) can be chosen as s{zo) = Ccl- 
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The reasons why we write a Umit in (4.23) will be seen in the proof of Theorem 2.3. 
At this point let us just say that we need to use Lemma 4.4 for the maximal choice 
s{zo) = Cei, in the cases when we know that ]D)p^^(zo) C ^ but do not know the same 
about the closure of D^^^ (^o). In light of continuity of z i-^ \ f{z)\, once s{zo) < Ccl, 
the limit is totally superfluous. 

Now we proceed to state a corresponding upper bound on 15(2) |: 

Lemma 4.5. There exists a constant A3 € (0, 00) and, for each C e (0, 00) and any 
sequence 7l obeying the assumptions (2.19), there exists a number L3 < 00 such that 

sup \g{z)\ < A3dL{zo)L''C{zo)''' (4.24) 

z: \z—zo\<C5l{zo) 

holds for any L> L3 and any zq S with ^c5l{zo){^o) C ^. 

With Lemmas 4.4^.5 in hand, the proof of Theorem 2.3 is rather straightforward. 

Proof of Theorem 2.3. Let m and n be distinct indices from TZ and let us abbreviate 
■^^j^ — '^y^{{m,n}) and = S^^{{m,n}). Let f{z) and g{z) be the functions 
from (4.20-4.21). Let Bi, B2, Ci, £2, C2 and A3 be the constants whose existence is 
guaranteed by Lemmas 4.3-4.5 and let Li be as in Lemma 4.3. Since A3 appears on the 
right-hand side of an upper bound, without loss of generality we can assume that 

C2A3 > Ci. (4.25) 

Further, let us choose the constants C and D such that 

C = Ci+C2A3, and D = Bi+B2. (4.26) 

Next, let L2 be the constant for which Lemma 4.4 holds for both C = C2 and C = 
C/A3 and for both = ^36"^^ and = AsL'^e-^'*^^^''. Finally, let L3 be the 
constant for which Lemma 4.5 holds with C as defined above. 

The statement of Theorem 2.3 involves two additional constants chosen as follows: 

First, a constant B for which we pick a number from (0, "^^1) (e.g, Bi/3 will do). 

Second, a constant Lq which we choose such that Lq > max{Li, L2, -^3} and that the 
bounds 

(4.27) 

hold true for all i > Lq. Fix L > Lq and consider the set 

^ = {^0 e %^ : Dc5^(.o)(zo) C (4.28) 

Notice that our choice of Lq guarantees that ^ C C ^k/(4l) n (.J^m U^„), while 

the fact that C < for both choices of C above ensures that for any zo £ the 

disc (zq) is contained in ff. These observations verify the assumptions of 

Lemma 4.4 — with cl = A36l{zo) and C equal to either C2 or C/A3 — as well as of 
Lemma 4.5, for any zq G 

First, we will attend to the proof of claim (2). Let zq G i?^ n ^ be a root of 
■^l' = 1^9- Lemma 4.4 with C = C2 and = ^315^(^0) and Lemma 4.5 then imply 
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the existence of a radius s{zo) with s{zo) < C2eL = C2A^5l{zo) < C5l{zq) such 
that 

|/(z)| > |.g(z)|, zeae,(zo) (4.29) 

holds for s — s(zo). (Note that here the limit in (4.23) can be omitted.) Hence, by 
Rouche's Theorem, / and f + g have an equal number of roots in I}>s(zo){zq), in- 
cluding multiplicity. In particular, the function / has a root zi in Ds(2„)(zo) which by 
Lemma 4.4 lies also in ^k/(2L)- Since s{zq) + Cie~'^^ < C(5l(zo) by the definition 
of C and the second bound in (4.27), we may use Lemma 4.3(3) to infer that the equa- 
tions (2.17-2.18) have a solution z E (zi) C I!>csL(zo)i^o)- Moreover, (4.27) 
implies that CSl{zo) < BiL^'^ so by Lemma 4.3(1) there is only one such solution in 
the entire disc ID)c5i(zo) (^o)- 

Next, we will prove claim (3). Let zq S i2L[Q) n be a solution to the equations 
(2.17-2.18). By Lemma 4.3(2), there exists a root zi e (zo) C I])csL(zn){za) 

of the function /. Lemma 4.3(1) then shows that zi is in fact the only root of / in 
'OcSl{zo){^o)- Applying Lemma 4.4 for the point zq and the choices cl — Az5l{zq) 
and C = CjA-i in conjunction with Lemma 4.5, there exists a radius s(zo) such that 
(4.29) holds true for any s < s(zq) sufficiently near s(zo). Moreover, by the bound 
(4.25) we know that z\ e ©^^^.^^(zo) C Dc2ei,(-2o) is a root of / within distance 
ci^L from zo, and so the last clause of Lemma 4.4 allows us to choose s(zo) = 
C(5l(zo). Let So < .s(zo) be such that (4.29) holds for s G (so, s(zo)) and pick an 
s G (so,s(zo)). Rouche's Theoremfor the discsDs(zo) and the fact that / has only one 
root in Dc5^(2(,)(zo) imply the existence of a unique zero z of /(z) + ^(z) = Z'^^{z) 
in Ds(zo). The proof is finished by taking the limit s \ C6l(zq). 

Further, we will pass to claim (4). Let z\ and Z2 be two distinct roots of Z^^^ in "^^^ 
such that both D^^-d(zi) C 6 and D^^-d(z2) C & are satisfied. We will suppose 
that \zi — Z2\ < BL~'^ and derive a contradiction. Let z — ^{zi + Z2) he the middle 
point of the segment between zi and Z2. Since |zi — Z2I < BL~'^, a simple geometrical 
argument shows that the disc of radius s = ^BL~'^ centered at z is entirely contained 
in ViBL-'^izi) U ^BL-''{^2) C ff. Next, by Lemmas 4.4-4.5, there exist two roots z[ 
and Z2 of / such that z[ e ^cs{zi){^i) aiid Z2 e ^C6{z2)i^2)- (We may have that 
zi = Z2, in which case Zi = Z2 would be a degenerate root of /.) Now our assumptions 
on B and Lq imply that 

^BL-''>^L-'' + CSl{zi) >\z-z,\ + \z,-z[\>\z-z[\, (4.30) 

and similarly for Zg. Consequently, both z'l and Z2 lie in Ds(z). But this contradicts 

Lemma 4.3 and the bound ^B < Bi, implying that ©^(zo) contains at most one 
non-degenerate root of /. 

Finally, we will prove claim (1). Let zq € ^ n %l{Q) with I])uL-d{z) C ^. 
According to Lemma 4.3(4), the disc D^^L '' (^) contains at least one one solution zi 
of the equations (2.17-2.18). Checking that ^2^"'' + CSl{zi) < {B2 + Bi)L-'^ in 
view of (4.27) and the definition of B, we know that Dp^^^^) (zi) C and we can use 
already proven claim (3) to get the existence of a root of Z'^^" in I^cSl^zi) {zi), which is 
a subset of D£)^-<i(zo). □ 

This concludes the proof of Theorem 2.3 subject to the viiUdity of Lemmas 4.3-4.5. 
The proofs of these lemmas have been deferred to Sect. 5.2. 
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4.3. Proof of Proposition 2.4. Fix distinct indices m,n ^ TZ. Our strategy is to first 
prove the claim for the density of the solutions of the equations (2.17-2.18), 

~p''^:n\^) = ^a\!^L{{m,n})f^^i,{z)\, (4.31) 

and then to argue that the density pm'n yields the same limit. 

Let Zo G '^{{m, n})\^, where ^ is the set of all multiple points. By Theorem 2.1 
and Assumptions A1-A2, there exists an e > such that, throughout the disc D,, = 
De(2;o) C 1^, we have Q(z) C {to, n} and the the function FTO,n (2) = Cm(^)/Cn(-z) is 
twice continuously differentiable and nonvanishing. Clearly, all solutions of the equa- 
tions (2.17-2.18) in De must lie in the set 

= e De: \Fra,n{z)\ = {(In / qmf ] ■ (4.32) 

Denoting the set ({m, njnlDie by '^'•■°°\ we now claim that for sufficiently small e, the 
sets can be viewed as differentiable parametric curves 7: f+) — > 

and 7^^) : (tL^^ , 4^') ^ D, for which 

(1) tL^' ^ t- andi^^' ^ t+ 

(2) 7^^' 7 uniformly on S 

(3) Vi — > V uniformly on [t- , t+) 

hold true as L ^ 00. Here VL{t) = ^7^^H*) ^nd y{t) = ^j{t) denote the tan- 
gent vectors to 7^^^ and 7, respectively. 

We will construct both curves as solutions to the differential equation 

dz{t) . d,<j) {z{t)) 

"dT ~ \d,(f>m,n{4m ' 

with (pm,n{z) = log|-FTO,„(2;)| (note that for e small enough, the right hand side is a 
well defined, continuously differentiable function of z{t) G ©c by Assumptions Al- 
A2 and the fact that \dz4'm,n{zo)\ > ct/^ according to Assumption A3). In order to 
define the curves 7^^^(-) and 7(-) we will choose a suitable starting point at t = 0. For 
7(-), this will just be the point zg, while for 7^^^(-) we will choose a point Zq^^ S 
which obeys the conditions (j)m.n{z'i^^) = riL and \zq — z'q^\ < Sa^^r/L, where 77^ = 

L~'^\og{qn/qrn)- To coustruct the point z^^^ G D^, we use again the smoothness of 
4>m,n- Namely, by Assumption Al-2, the function (f)m.n{x + iy) = log \Fm,n{x + iy)\ is 
twice continuously differentiable on if e is sufficiently small, and by Assumption A3 
we either have \d4>m,n{x + w)/9x\ > a/3, or |9(^„_„(a; + iy)/dy\ > a/3. Assuming, 
without loss of generality, that \d(i>rn,nix + iy) / dy\ > a/3 on all of De, we then define 
Zq^'' to be the unique point for which Rczq^-* = Rezo and (f>m.n{zQ^'') = 77^. By the 
assumption |i9</)TO^„(a;+iy)/52/| > a/3, wethenhave |zo — -^o^^l < 3a~"^r7L, as desired. 

Having chosen Zq"\ the desired curves j^^^ : (t^^ , ) —^ and 7 : {t-,t+) —^ 
Pe are obtained as the solutions of the equation (4.33) with initial condition 7(^^(0) = 
Zq"^ and 7(0) = zq, respectively. Here t^^\ t-, and t+ are determined by the 
condition that t^:^^ and t_ are the largest values t < for which j^'"\t) G 9De and 
7(t) G 3De, respectively, and t[^^ and t+ are the smallest values t > for which 
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J^^\t) G 9De and 7(f) e 9Be, respectively. Since the right-hand side of (4.33) has 
modulus one, both curves are parametrized by the arc-length. Moreover, decreasing e 
if necessary, the functions 7^^' can be extended to all t € To see that the 

limits in (1-3) above hold, we just refer to the Lipschitz continuity of the right hand 
side of (4.33) and the fact that, by definition, |7'^H0) -7(0)| = 0(L-'').Letirbethe 
Lipschitz constant of the right-hand side of (4.33) in a neighborhood containing 7^^) {t) 
for all t e (t_ , Choosing e so small that both t+ — t_ and t^^^^ — f*^-* are less than, 
say, 1/{2K), integrating (4.33) and invoking the Lipschitz continuity, we get 

sup M^Ht)-7WI < l7^^nO)-7(0)| + | sup |7(^)(i)-7(i)l- (4-34) 

t-<t<t+ t-<t<t+ 

This shows that 7^^^(f) 7(t) uniformly in f e {t-,t+). Using Lipschitz continuity 
once more, we get a similar bound on the derivatives. But then also the arc-lengths 
corresponding to 7^^^ must converge to the arc-length of 7, which shows that also 

Consider now the curve ^{t) . Given that | Fm,n {z)\is constant along 7, we have 

i^i^Ml) = ilisM)) = ,„gF„,.(.)U,,„,(,). (4.35, 

Referring to Assumption A3 and the fact that \v{t)\ = 1, we find that the modulus of the 
left-hand side is bounded below by a. Using continuity of the derivative ^ ArgF^ „ 

in Dg, we observe that one of the two alternatives occurs on all the interval 4^^): 

dArgJ-^,„(7W(^)) a d Argf^,„(7W(^)) a 

either ^-^ > — or ^— ^ < , 

dt - 2 dt - 2' 

(4.36) 

provided e is sufficiently small. By Lemma 4.3, the disc contains a finite number 
k = 2eL''-pm'n (zq) of solutions of the equations (2. 17) and (2. 18) which in the present 
notation read 

\Fm,n{z)\ = (^) , (4.37) 

L'^AvgFm,n{z) = 7rmod27r. (4.38) 

Assuming, without loss of generality, that the former alternative in (4.36) takes place, 
and ordering all the solutions consecutively along the curve 7*^^\ i.e., letting Zi = 

7(^)(ti),..., ^fe = 7(^)(ife) where iL^^ <ti< ■■■ <tk <t^^\wehave 

ArgF„,„(^,+i) - ArgF„,„(^,) = 2nL-'^ (4.39) 
for any j = 1, . . . , fc — 1, as well as 

AigFrr,,n{zi) " ArgJ^„,„(2_) < 2nL-'^ (4.40) 

and 

Arg^;„,„(^+) - Aig Frr,,n{zk) < 2ttL-'^. (4.41) 
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In view of the first equality in (4.35) rephrased for 7^^\ the left hand side of (4.39) can 
be rewritten as 



ArgF„,„(2j+i) - ArgFm,n{zj) = / 
and thus 



dlogF„,„(7(^)(t)) 



+ 



dlog^;„,„(7(^)(0) 



dt 



dt - 2fc7rL-<* 



dt 



< 2-kL- 



dt 



(4.42) 



(4.43) 



Let us divide the whole expression by L''- and take the limit i ^ oo. Now 7'^^ converge 
to 7 along with their first derivatives, uniformly in t e {t- , t+), and the limits fjf' con- 
verge to t± . The Bounded Convergence Theorem then shows that the integral in (4.43) 

converges to a corresponding integral over 7. Recalling that pm,'n (zo) = k/ {2eL'^), we 
thus get 



L— >oo ' 47re 



dlogF„,„(7o(t)) 



dt 



dt 



= -r— \dAog Fm,n{z)\\dz\ 
47re 



(4.44) 



where the last integral denotes the integration with respect to the arc length. Taking into 

account the Lipschitz continuity of \dz \ogFm,n{z)\, the last integral in (4.44) can be 
approximated by [\dz \ogFm,n{zo) \ + 0{e)) |7|. By the smoothness of the curve 7, we 
estimate its length by I7I = 2e(l + 0(e)), so that 



^}J^}^^^r^,ni^0) = ^\9zlogFm,n{zo)\ = ^ 



2w 



dzCmizo) dzCnizo) 



Cmizo) Cn{zo) 



(4.45) 



To finish the proof, we need to show that pm^'n (zo) will converge to the same limit. 
According to Theorem 2.3, we have 

\\nin'D^{z)\ - \ f2L{{m,n})nI]>^{z)\ \ < 2 (4.46) 

for all ^ G $f(m, n) such that | Q{z)\ = 2 and e sufficiently small. Hence 

\pl^;n\z)-pl^,i\z)\<^, (4.47) 
and the claim of the proposition follows by (4.45). □ 



4.4. Multiple phase coexistence. In this section we wiU prove Theorem 2.5, which 
deals with the zeros of Z^'^ in the vicinity of multiple points. Let zm E ^ be a multiple 
point and let Q = Q{zm)- For each m e Q, let 4>m{L) and be as in (2.23). Define 
the functions 



~g{z) = Zr(^)C(^M)-^' - f{z), 



(4.49) 
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and 

^{z) = exp{maLxRe{vm{z - zm))}- (4.50) 

mGQ 

As in the case of two-phase coexistence, the proof uses Rouche's Theorem for the 
functions / and f + g. For this we will need a lower bound on |/| and an upper bound 
on |^|. 

Lemma 4.6. Suppose Assumptions A and B hold. Given Q C TZ with \Q\ > 3 and 
abbreviating q = \ Q\ and Rl = let (cl) be a sequence of positive numbers 

such that 

lim L'^'^CL = oo but lim L'^'^-'^/Icl = 0. (4.51) 

Then there is a constant L5 < oo such that for any zo & C and any L > L5 there 
exists s{zo) S [Rl/QjRl] for which the bound 

inf |/(z)| >L'^eiC(^of' (4.52) 

z: \z—zo\=s(zo) 

holds. 

Lemma 4.7. Let zm G & be a multiple point, let Q = Q{zm), q = \Q\, and Rl = 
jiigj-g exists a constant e (0, 00) and, for each sequence (pl) of posi- 
tive numbers obeying (2.25), a number Lq < 00 such that if L > Lq then P^'^ (-2m) C 
'^k/l{Q)' where p'^= pl + Rl- Furthermore, we have 

sup \g{z)\<AQplL%zo)''' (4.53) 

z : \z—zo\<Rl 

whenever zq G ]S>p^ {zm)- 

With these two lemmas we can proceed directly to the proof of Theorem 2.5. 

Proof of Theorem 2.5. The proof is close in spirit to the proof of Theorem 2.3. Let zm 
be a multiple point and let Q = Q{zm). Consider a sequence (pl) of positive num- 
bers such that (2.25) holds. Choosing = ^epi. where Aq is the constant from 
Lemma 4.7, we note that the conditions (4.51) are satisfied due to our conditions on pL 
from (2.25). We will then prove Theorem 2.5 with Lq = max{_L5, Lq}, where L5 and 
Lq are the constants from Lemma 4.6 and 4.7, respectively. The proof again boils down 
to a straightforward application of Rouche's Theorem. 

Indeed, let L > Lq and note that by Lemmas 4.6 and 4.7, for each zq € ©pj^ (zm) 
there is an s{zo) e [RL/q, Rl] such that on Ds(2j,)(zo), we have 

\fiz)\ > \g{z)\. (4.54) 

Consider the set of these discs Vig(^zo){^o) — one for every zq € Dp^ (^^m). These discs 
cover the closure of Dp^ (zm), so we can choose a finite subcover S. Next we note that 
(4.54) implies that neither / nor f + g have more than finitely many zeros in EDp^ (zm) 
(otherwise, one of these functions would be identically zero). Without loss of generality, 
we can thus assume that the discs centered at the zeros of / and f + g in Dp^^ (zm) are 
included in §. Defining = UdgS clearly have Bp^- (zm) C C Dp' (^m). 

Let now 3C be the set of all components of \ IJugg dH). Let J(f E 3C he one 
such component. By (4.54) we know that \f{z)\ > \g{z) \ on the boundary of and 
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Rouche's Theorem then guarantees that / has as many zeros in as f + g, provided 
we count multiplicity correctly. Moreover, both functions / or / + g have no zeros on 
Upgs^B. Since = '(z)C(^m)-^" andC(zM)-^' > 0, the zeros of /+ 5 
are exactly those of The above construction of ^ and S then directly implies the 
desired correspondence of the zeros. Namely, in each G 3C, both / and Z^'^'^ have 
the same (finite) number of zeros, which can therefore be assigned to each other. Now 
/ and have no zeros in ^ \ Ujtgsc choosing one such assignment in each 

J?r G 3C extends into a one-to-one assignment of f2j^ n ^ and fiiiQ) H Moreover, 
if z e H'ln J(f and z G fl for some ^ £ X (which is required if z and 

z are the corresponding roots), then z belongs to the disc D e S centered at z and z 
belongs to the disc D e S centered at z. Consequently, z and z are not farther apart than 
Rl = L-''(i+i/«). This completes the proof. □ 



4.5. Proof of Proposition 2.6. Assuming that L lol < 7l, it clearly suffices to ascer- 
tain that 

y yj^{Q)n^C y B>p,{zm). (4.55) 
C: |Q|>3 zuesn.^ 

To this end let us first observe that continuity of the functions implies 

lim ^^JQ) = fl ^„ (4.56) 

i— >oo ' ' 

meQ 

since 7/, 0. The set ^ fl ^ is finite according to Theorem 2.1. Hence, there exists 
a constant (5o > and, for each S e (0, 60], a constant Lq = Lo{S), such that the discs 
^^(^m), -Zm G ^ n are mutually disjoint, 

Q{z) C Q{zm) whenever z e Bs{zm), (4.57) 

and 

y ^7JQ)n^C y ©5(0m) (4.58) 

Q: |C|>3 ZMG^n^ 

whenever < S < So and L > Lo(<5). It is therefore enough to show that there exist 
constants x > ^nd 5 G (0, (5o) such that for any multiple point zm G ^, we have 

Bs{zm) n y^AQizM)) C Bp^(^M) (4.59) 

once PL > xiL and L > Lo{5). 

We will prove (4.59) in two steps: First we will show that there is a constant x > 
such that for any multiple point zm, any z ^ zm, and any n € Q{zu), there exists 
m e Qizyi) for which 

Re[(z - ZM){vn{zM) - Vynizu))] > 2xk - | , (4.60) 

and then we will show that (4.60) implies (4.59). To prove (4.60), we first refer to the 
fact that we are dealing with a finite number of strictly convex polygons with vertices 
{vkizu) '■ k e Q{zm)} according to Assumption A4 and thus, given z and n, the label 
m can be always chosen so that the angle between the complex numbers z — zm and 
Vnizu) — u„, (zm) is not smaller than a given fixed value. Combining this fact with the 
lower bound from Assumption A3, we get (4.60). 
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We are left with the proof of (4.59). Let us thus consider a multiple point zm G ^ 
with Q{zm) = Q, and a point z G (zm) \ Dp^ (zm). We wiU have to show that there 
exists an m e Q with z ^ .9'^^ (m). Recalhng that Q(z') C Q for all z' e ©^(zm), 
let n G Q be such that iCn (2) I = CC-^)- Choosing m G Q,{zyy) so that (4.60) is satisfied 
and using, as in the proof of Lemma 3.1, Fn^^iz) to denote the function = 
Cn(2^)/CTO(-2). we apply, as in (3.9), the Taylor expansion to log |-Fn,m(-2)| to get 

log \Fn,m{z)\ = Re[(2; - Zyi){Vn{zyi) - Vm{zyi)f\ + 0(\z ~ ZmP) 

>xV-zm\>XPl- (4.61) 

Here we also used that \Fn^m{zM)\ = 1 and assumed that b was chosen small enough 
to guarantee that the error term is smaller than x| -2 — -^m | • As a result, we get 

|Cm(^)| < e-^P^^{z) < e-^^C(^) (4.62) 

implying that z ^ (m). Thus, the inclusion (4.59) is verified and (4.55) follows. 
□ 



5. Technical lemmas 



The goal of this section is to provide the proofs of Lemmas 4.3-4.7. We will begin with 
some preparatory statements concerning Lipschitz continuity of the Cm and C,. 



5.1. Lipschitz properties of the functions log |Cm| and logC- In this section, we prove 
two auxiliary lemmas needed for the proofs of our main theorems. For any 2:1, 02 G C, 
we will use \z\ , Z2] to denote the closed segment 

[zi,Z2] = {tei + (l-t)z2: tG [0,1]}. (5.1) 

The following Lipschitz bounds are (more or less) a direct consequence of formulas 
(2.9) and (2. 1 1) in Assumption B. 

Lemma 5.1. Suppose Assumptions A and B hold and let k, t, and M be as in Assump- 
tion B. Let m G TZ, and let zi,Z2 G ^k/l("^) be such that [-21,-22] C '^k/l{it^)- 
Then 



Cm(zi) 



Cm (^2) 

Moreover, for all zy^ziG G such that \z\., 02] C G, we have 

C(2:i) ^ „M|zi-;j2| 



C(^2) 



< e'" 1^1-^=1. (5.3) 



Remark 7. Since z ^ \C,ra{z)\ are all twice continuously differentiable and hence Lip- 
schitz throughout G , so is their maximum z ^ C,(z). The reason why we provide a 
(rather demanding) proof of (5.3) is that we need this bound to hold uniformly through- 
out G and the constant M from Assumption B(3) to appear expUcitly on the right-hand 
side. The first part of the lemma underlines what is hard about the second part: On the 
basis of Assumption B, the uniform Lipchitz bound in (5.2) can be guaranteed only in 
the region where m is "almost stable." 
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Proof of Lemma 5.1. Let [zi, Z2] C ^^/^(m). The bound (5.2) is directly proved by 
combining (2.9) with the estimate 



|l0g|C(f)(^l)| - l0g\C^^\z2)\\ < M\Z, - Z2\, 

impUedby (2.11). Indeed, introducing (^(t) = Cm\zi + t{z2 — zi)), we have 



1 d\ip{t)\ 



if{t) At 

implying (5.4). By passing to the limit L 



< 



1 



dt 



<M\Z2-Zi\ 



00, we conclude that 



|l0gC(2:i) - l0gC(-22)| < M\zi - Z2\ 



(5.4) 



(5.5) 



(5.6) 



holds provided [z\,Z2\ C S^m- 

To prove (5.3), let zi, Z2 G with [zi, Z2] C l^. If the segment [zi, Z2] intersects the 
coexistence set ^ only in a finite number of points, then (5.3) is an easy consequence 
of (5.6). However, this may not always be the case and hence we need a more general 
argument. Note that continuity of both sides requires us to prove (5.3) only for a dense 
set of points zi and Z2. This and the fact that each compact subset of contains only a 
finite number of multiple points from ^ = {z € l^: \ Q{z) \ > 3} permit us to assume 
that zi, Z2 ^ and that the segment [zi, Z2] does not contain a multiple point, i.e., 

[zi,z2] = 0. 

Suppose now that the bound (5.3) fails. We claim that then there exist a point x G 
[zi, Z2], with X ^ Zi, Z2, and two sequences and of points from [zi,x] (1 
and [a;, Z2] fl 'S, respectively, such that the following holds: 

(1) Xn 7^ Vn for all n and lim„^oo a^n = lim„^<x) Vn = x. 

(2) There exists a number M' > M such that 



avn) 



> M'\x„ - y„| 



(5.7) 



for all n. 



The proof of these facts will be simplified by introducing the Lipschitz ratio, which for 
any pair of distinct numbers x,y [z\, Z2] is defined by the formula 



R{x,y) 



\x-y\ 



(5.8) 



The significance of this quantity stems from its behavior under subdivisions of the in- 
terval. Namely, if x and y are distinct points and z <£ (x, y), then we have 



R{x, y) < inax{R{x, z),R{z, y)}, 



(5.9) 



with the inequality being strict unless R{x, z) = i?(z, y). 

To prove the existence of sequences satisfying (1) and (2) above, we need a few 
observations: First, we note that M' = R{zi, Z2) > M from our assumption that (5.3) 
fails. Second, whenever x, y e [zi, Z2] are such that R{x, y) > M, then (5.6) implies 
the existence of x',y' e [x, y] such that x', y' and R{x' , y') > R{x, y). Indeed, we 
choose x' to be the nearest point to x from the closed set [x, y] n ^, and similarly for y' . 
The fact that the Lipschitz ratio increases in the process is a direct consequence of (5.9). 
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Finally, if distinct x, y G [zi, Z2] H satisfy R{x, y) > M, then there exists a pair of 
distinct points G [a;,y]n^^ such that [x' — y'l < ^\x — y\ and R{x' ,y') > R{x,y). 
To prove this we use (5.9) with z = ^{x + y) to choose the one of the segments [x, z] 
or [z, y] that has the Lipschitz ratio not smaller than R{x, y) and then use the preceding 
observation on the chosen segment. 

Equipped with these observations, we are ready to prove the existence of the desired 
sequences. Starting with the second observation above applied for x = z\ and y =- 
Z2, we get xi,X2 G [z\,Z2] n such that R{x\,X2) > M' . Notice that xi ^ zi 
and X2 =/= Z2 since z\,Z2 ^ ^. Next, whenever the pair x„, ?/„ is chosen, we use the 
third observation to construct the pair a:;„+i, y„+i G \xm y^] n 5f of points such that 
|a;„+i - y„+i| < \\xn - J/nl and ii!(a;„+i, > R{xn,yn) > M' . Clearly, the 
sequences (a;„) and (y„) converge to a conmion limit x G [a;i,t/i], which is distinct 
from Zx and Z2. 

We will now show that (5.7) still leads to a contradiction with (5.3). First we note 
that the point x, being a limit of points from \ ./#, is a two-phase coexistence point 
and so Theorem 2.1(2) applies in a disc 0^{x) for e > sufficiently small. Hence, 
there is a unique smooth coexistence curve connecting x to the boundary of De(x) 
and, since (x„) and (y„) eventually lie on its tangent vector at x is colinear with the 
segment [zi, 2:2]. Since in (.x), the coexistence curve is at least twice continuously 
differentiable, the tangent vector to '6' has a bounded derivative throughout IDf.(.T). As a 
consequence, in the disc (x) with 5 <e, the curve will not divert from the segment 
[01, 02] by more than CS"^, where C = C(e) < 00. 

Now we are ready to derive the anticipated contradiction: Fix n and let 5n be the 
maximum of — a;| and |y„ — x\. Let e be a unit vector orthogonal to the segment 
[zi, Z2] and consider the shifted points x'^ = Xn + ^CS^e and y'n= yn + 2C^^e. Then 
we can write 

^ C(£n) CK) C(K) ..10^ 

Civn) C«)C(y;)C(2/n)' ^ ■ ' 



Assuming that n is sufficiently large to ensure that 5„ ^/l + 4CPSf^ < e, the segment 
[x'^jy'^] lies in 3^{x) entirely on one "side" of and is thus contained in J^^n for 
some m G 7?.. On the other hand, given the bounded derivative of the tangent vector 
to each segment [xn,x'^] and [t/„, y'^] intersects the curve ^ exactly once, which in 
hght of Xn,yn G $f happens at the endpoint. This means that also [xn, x'^] C S^m and 
[Vn, y'n] C for the same m. Consequently, all three ratios can be estimated using 
(5.3), yielding 

R{xn,yn) < M ~ + ~ + ~ < M + AMC6n, (5.11) 

where we used that \x'^ — y'^\ = |x„ — y„| and |x„ — y„| > Sn- But 5„ — > with 
n ^ oo and thus the ratio R{xn, yn) is eventually strictly less than M' , in contradiction 
with (5.7). Hence, (5.3) must have been true after all. □ 

The previous lemma will be particularly useful in terms of the following corollary. 

Corollary 5.2. Suppose that Assumptions A and B hold and letO < k < k, where k is 

the constant from Assumption B. Then there exist constants c < oo and L4 < oo such 
that the following is true for all L > L4 and all s < c/ L: 

(1) Form G TZ and z G >^s/(2L)(™) with ©^(z) c ^, we have 

B,{z) C ^s/lM. (5.12) 
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(2) Forz€ ^ with Iis{z) C 0, the set 

Q' = {meTe: C^«/L(m)} (5.13) 

in non-empty and 

D,,(z) C f/s/L(Q')- (5.14) 

(3) For 71, < k;/(2L), Q c TZ and z € %l{Q) ^ '^2k/i,(Q) with Bs{z) C ^, we 
have 

^s{z) C %i(Q). (5.15) 

Proof. Let M be as in Assumption B. We then choose c > sufficiently small and 
1^4 < oo sufficiently large to ensure that for L > L4 we have 

TTTT - -TjLe-''^ > 2c. (5.16) 
8M M - 

First, we will show that the claims (1), (2), and (3) above reduce to the following 
statement valid for each m G 71: If z, z' are complex numbers such that the bound 
\z — z'\ < 2c/L, the inclusion [z, z') C ff, and z ^ G\ >9^K/Lijn) hold, then also 



[z,z')^&\y-^,i2L){m). (5.17) 

We proceed with the proof of (1-3) given this claim; the inclusion (5.17) will be estab- 
lished at the end of this proof. 

Ad(l): Let z e -^7(2L) with Dig{z) C and assume that (5.12) fails. Then there exist 
some z' e i^' \ /L (™) with \z — z'\ < s and [z, z'\ C &. But by (5.17), this implies 

[z', z)n^£/(2L)(TO) = 0, which means that [z', z]n^;5/(2i,)(m) = 0. This contradicts 
the fact that z G ■5^K,l{iL)ijn). 

Ad (2): Let z e ^ with ©^(z) C (ff. By the definition of stable phases, there is at 
least one m E TZ such that z e C ^ji/(2L){fn)- Combined with (5.12), this 
proves that the set Q' is non-empty. To prove (5.14), it remains to show that Ds(z) C 
^ \ ^ki(iL) (to) whenever m ^ Q'. By the definition of Q', m ^ Q' impUes that there 
exists a z' G ID)s(z) such that z' G \ S^a/hii^)- Consider an arbitrary z" G 91D)s(z). 
For such a z", we have that |z' — z"| < 2c/ L and [z',z") C so by (5.17), we 
conclude that [z', z") C ^\=5^s/(2L)(m). Since this is true for all z" e 9Ds(z), weget 
the desired statement Ps(2) C ^ \ yk/(2L) (w). 

Ad (3): Let Q C TZ, z e %^iQ) n '^2ft/L(S) and D,(z) C ^. If m € S, then 
z € .y-yj^ (m) C ^K/(2L)('Ti) by the definition of '^^^ (Q) and the condition that 7/, < 
k/{2L). With the help of (5.12), this implies that P^z) C YR/Lim) for all m e Q. 
Recalling the definition of '^^/^(Q), we are left with the proof that ©^(z) d \ 

■^K/{2L){in) whenever m ^ Q. But if m ^ Q, then z & \ S^a/Li'm) because we 
assumed that z e '^2r/l{Q)- By (5.17) we conclude that [z, z') C ^ \ ^k/(2l)(to) 
whenever z' e 9Ds(z), which proves ©^(z) C ^ \ '5^k/(2l)(to)- 

We are left with the proof of (5.17), which will be done by contradiction. Assume 
thus that m G TZ and let z, z' be two points such that |z — z'| < 2c/ L, [z, z') C ^ and 
z e ^ \ ^s/l(™) hold, while (5.17) fails to hold, so that [z, z') n J=^s/(2Z,) (m) 0- 
Letzi G [z, z')n,yj^/(2L)(m). Since [z.z') C i^?, we have in particular that [zi , z] C ^■ 
Let Z2 be defined as the nearest point to zi on the linear segment [zi,z] such that 
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zi ^ -5l3ft/(4L)(™)- By continuity of the functions (fe, we have [^1,22] C ^k/l(™) C 
'5^K/L{m) so that the bounds in Lemma 5.1 are at our disposal. Putting (5.2-5.3) to- 
gether, we have 



Cm(^l) C(^2) 



C(^l) 



C,m{z2) 



< e 



2e-'^^+2M|zi-22| 



(5.18) 



Now, since z\ € '^k/{2L){fn) and ^ ■^3k/(4l)('7i)> we can infer that the left-hand 
side is larger than e**/ ^^^^ . Hence, we must have 



1-21 - -221 > 



At _ 1 . 

8MZ ~ 



-tL 



> 



2c 



(5.19) 



where the last inequality is a consequence of (5.16). Now zi, Z2 G [z, z') implies \zi — 
Z2\ < \z — z'\, which contradicts the assumption that \z — z'\ < 2c/ L and thus proves 
(5.17). □ 



5.2. Proofs of Lemmas 4.3-4.5. Here we will establish the three technical lemmas on 
which the proof of Theorem 2.3 was based. Throughout this section, we fix distinct 
m,n & TZ and introduce the abbreviations = ■^e{{'m', n}) and = '^eiim, n}). 
We will also let / and g be the functions defined in (4.20^.21). 

First we will need to establish a few standard facts concerning the local inversion of 
analytic maps and its behavior under perturbations by continuous functions. The proof 
is based on Brouwer's Fixed Point Theorem, see e.g. [30, Chapter 2]. 

Lemma 5.3. Let zq ^ C, e > 0, and let (f) : Of (zo) ^ C be an analytic map for which 

\cl>'{zo)\-'\<l>\z)-<l>'{zo)\<\ (5.20) 

holds for all z G ©^(zo). Let 5 < €\(j)'{zo)\/2. Then, for every w G ©^(^(zo)). there 
exists a unique point z E (zo) such that (biz) = w. 

In addition, let rj G [0, 5/2) and let 6 : Dt(zo) C be a continuous map satisfying 

\0{z)\<r^, zeD,(zo). (5.21) 

Then for each z € De(zo) with <p{z) € V>fj{(j){zo)) there exists a point z' € De(.2o) such 
that 

cj>{z')+e{z')=cf>{z). (5.22) 
Moreover, \z' — z\ < 2ri\(p' {zo)\~^ ■ 

Proof. Following standard proofs of the theorem about local inversion of differentiable 
maps (see, e.g., [13], Sect. 3.1.1), we seek the inverse of w as a fixed point of the (ana- 
lytic) function z il){z) = z + (j)' {zq)"^ {w — 4>{z)). The condition (5.20) guarantees 
that z I— > tp{z) is a contraction on Iii^{zo). Indeed, for every z G ^e{zo) we have 

|V'(z)| = |1 - ct>'{zo)-'ct>'{z)\ < \^'{zo)\-'\ct>'{z) - ^'{zo)\ < i, (5.23) 

which implies that |V'(z) — V'(^') I ^ ^\z — z'\ior all z, z' G ©^(zo)- The actual solution 
to (p{z) = u> is obtained as the limit z = lim„_»oo z„ of iterations z„+i = ip{zn) 
starting at zq. In view of the above estimates, we have \zn+i — Zn\ < \\zn — z„_i| 
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and, summing over n, we get \zn — zo\ < 2\zi — zo\ < 2|0'(zo)|^^|w — (t){zo)\. Since 
|u; — (l){zo)\ < 5, we have that Zn as well as its hmit belongs to De(^o)- 

Next we shall attend to the second part of the claim. The above argument allows 
us to define the left inverse of (p as the function (p~^ : I])s{(f){zo)) De(zo) such that 
(j)~^{w) is the unique value z e De(zo) for which = w. Let r] G [0,6/2) and let 
z e I!>e{zo) be such that ?!>(^) e I])n{4>{zo)). Consider the function tZ^ : ©^(^(^o)) ^ C 
defined by 

If (w) = 4>{z) - 9{(j>-^{w)). (5.24) 

By our choice of z and (5.21), we have |!f'('u;)| < 2ri for any w G 'Ds{(j){zo)). Thus, "f" 
maps the closed disc Ii2n{^{zo)) into itself and, in light of continuity of ^, Brouwer's 
Theorem implies that ^ has a fixed point w' in I}2ri{(f>{zo))- From the relation ^{u/) = 
w' we then easily show that (5.22) holds for z' = 4)^^{w'). To control the distance be- 
tween z and z' , we just note that the above Lipschitz bound on ijj allows us to conclude 
that \z' - z\< 2\(j)'{zo)\-'^\(l){z') - (j){z)\. Applying (5.22) and (5.21), the right-hand 
side is bounded by 2?7|(?i(2:o)|~^. □ 

Now we are ready to start proving Lemmas 4.3-4.5. The first claim to prove concerns 
the relation of the solutions of (2.17-2.18) and the roots of the function / defined in 
(4.20). 

Proof of Lemma 4.3. Let a, M and r be the constants from Assumption B. Let c and 
be the constants from Corollary 5.2 with k = k. The proof will be carried out for the 
constants Bi, C\ and Li chosen as follows: We let 

B,= ' ^,^16 + 4|log(Wa>n)l ^^^10^ (5.25) 

4M a a 

and assume that Li is so large that Li > L4 and for all L > Li, we have Cic"^^ < 
BiL~'^ and the bounds: 

^B.+B,)L-'^<'-<^, 2e-^ + ^<i (5.26) 

^(M + M2)(Bi+B2)i-'^ < J, (5.27) 
a z 

and also 

26"^-^ + 2MB1L-'' < L-'^, a > 2V2e-''^, (5.28) 

ttL-'^ + 2e-^^ < AL-'^ and Cie'^^ < iSaL"''. (5.29) 

Let us fix a value L > Li and choose a point zq e .^k/(2l) and a number s < {Bi + 
B2)L~''' such that Bsiza) C ^. Corollary 5.2(1) combined with the first bound in 
(5.26) implies that ©^(zo) C S^^/l- 

We will apply Lemma 5.3 for suitable choices of (p and defined in terms of the 

functions Fm^n ■ ID's(^o) — * C and Fm ji ■ Ds(zo) — > C defined by 

FmA^)=^f^ and fW(.) = ^. (5.30) 

We will want to define (j){z) as the logarithm of Fm}n{z), and 6{z) as the logarithm of 
the ratio Fm}i {z) / Fm,n {z), but in order to do so, we will have to specify the branch of 
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the complex logarithm we are using. To this end, we will first analyze the image of the 
functions F^}n{z) and F^^Xiz) I F„i.n{z). 

According to Assumption B2, for any z e Bs(-Zo) C ^k/l, we have \Frn}i{z)\ G 
(2/3, 3/2) in view of the second bound in (5.26) with the observation that |; < log |. 
A simple calculation and the bound (2.11) show that AvgFin}n{z) and ArgFm}i{zo) 
differ by less than 2M(Bi + 82)1-'^ < \. Indeed, the difference Arg F^l(z) - 

ArgF^^(zo) is expressed in terms of the integral of d^Fm^i/ Fm}n along any path in 
Ds(zo) connecting zq and z. The latter logarithmic derivative is bounded uniformly by 

2M throughout ©^(zo). Consequently, z ^ Fm}i{z) maps ©^(zo) into the open set 
of complex numbers {pe"^ : p G (|, |), |cj — cjo| < ^}, where = ArgFm,„(zo). 

The function z Fm}n{z) / Fm,n{z), on the other hand, maps Ds(zo) into the open 
set of complex numbers {pe"^ :pS < ^j, as can be easily inferred from 

Assumption B2 and the second bound in (5.26). Given these observations, we choose 
the branch of the complex logarithm with cut along the ray {re^"^°/^ : r > 0}, and 
define 

ci>{z) =\ogF^l{z) (5.31) 

and 

F^"-^ {z^ 

e{z) = \og'-^^ (5.32) 

Having defined the functions (j) and 0, we note that, by Assumptions A and B, (j) is 
analytic while 9 is twice continuously differentiable throughout Hg^zo). Moreover, 
these functions are directly related to the equations /(z) = and (2.17-2.18). In- 
deed, /(z) = holds for some z e Bs(zo) if and only if Fm}i{z) is an L'^-th root 
of —{qn/qm), i-C-, 'f'iz) = (log(q'n/9m) + i7r(2A; + 1))L^'^ for some integer k. Simi- 
larly, z e Bs(zo) is a solution of (2.17-2.18) if and only if (p{z) + 9{z) is of the form 
{log{qn/qm) + i'JT{2k + 1))L~'^ for some integer k. Furthermore, these functions obey 
the bounds 

a<\(t>'iz)\<2M, |(/)'(z)-<^'(2o)| <2(M + M")(Bi+S2)L-^ (5.33) 

and 

\eiz)\ <2e-^^, \9{z)-9{z')\<2V2e-^^\z- z'\ (5.34) 

for all z, z' G Bs (zq). Here the first three bounds are obvious consequences of Assump- 
tion B, while the third follows from Assumption B by observing that the derivative ma- 
trix D9{z) is bounded in norm by 2v^ times the right hand side of (2.10). Note that, 
in Ught of (5.26-5.27), these bounds directly verify the assumptions (5.20) and (5.21) 
of Lemma 5.3 for i] = 2e~'^^ and any e < s. We proceed by applying Lemma 5.3 with 
different choices of e to give the proof of (2-4) of Lemma 4.3, while part (1) turns out 
to be a direct consequence of the bounds (5.33-5.34). 

Indeed, let us first show that for s < BiL~'^ the disc Bs(zo) contains at most one 
solution to (2.17-2.18) and at most one root of the equation /(z) — 0. We will prove 
both statements by contradiction. Starting with the solutions to (2.17-2.18), let us thus 
assume that Zi,Z2 G Dg(zo) are two distinct solutions to the equations (2.17-2.18). 
Setting w\ = ^(zi) -|- 9{zi) and W2 = 4'{z2) + 9{z2) this means that w\ — W2 is 
an integer multiple of 27riL~''. However, the bounds (5.33) and (5.34) combined with 
the first bound in (5.28) guarantee that \wi - W2I < 4e~'"'^ + AMBiL^'^ < 21-"^ 
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and thus wi ~ W2- But then the bound \(t>{zi) — (f'{z2)\ > o:\zi — Z2I implies that 
\0{zi) — ^(22)1 > ct\zi — Z2\, which, in view of the second bound in (5 .28), contradicts 
the second bound in (5.34). Hence, we must have had zi = Z2 in the first place. Turning 
to the equation f{z) — 0, let us now assume that zi and Z2 are two different roots of this 
equation. Setting wi = (j){zi) and W2 = 4'{z2), we again have wi = W2, this time by 
the first bound in (5.33) and the very definition of Bi, which implies that AMBi = 1. 
But once we have wi = W2, we must have zi = Z2 since |</'(zi) — (6(22) | ^ o/\zi — Z2\ 
by our lower bound on (j>'{z), implying that there exists at most one z G Ds(zo) that 
solves the equation f{z) = 0. If such a solution z exists, Assumption B inomediately 
imphes that f'{z) ^ 0, and so 2 is a non-degenerate root of /. 

Next, we will show that within a C'ie~'^^-neighborhood of each solution zq of the 
equations (2.17-2.18) there is a root of /. Indeed, let e = CiC"^^ and S = 5e~^^. By 
the first bound in (5.33) and our choice of Ci, we then have 6 < €\(j)'{zo)\/2, so the 
first part of Lemma 5.3 is at our disposal. Since zq is assumed to be a solution to (2.17- 
2.18), we have that (/)(zo) + (^{zo) is of the form {\og(qn/qm) + i7r(2fc + l))L^'^, where 
fc is an integer. In light of the bound |^(2o)| < 2e~'^^, the disc IS>s{4>{zo)) contains the 
point ti; = (/)(2:o) + ^(-2o)-By the first part ofLemma5. 3, there existsapoint2 G ^eizo) 
such that (t){z) ~ w, implying that z is a root of /. 

As a third step we will prove that if zq is a root of /, then there exists a solution 
to (2.17-2.18) in D^^^-^i, (^;o). By the relation between / and <p we now know that 
(j){zo) is of the form {\og{qn/qm) + i7r(2A: + 1))L~'^ for some integer k. We again set 
e = Cie~'^^ and S = 5e~^^. Choosing t] = 2e~^^ and noting that 2r] < 6, we apply 
the second part of Lemma 5.3 to conclude that there must be a point z' G Oe{zo) such 
that (j>{z') + 0{z') = <l>{zo) = {log{qn/qm) + i7r(2fc + 1))L-^, which means that z' is 
a solution to (2.17-2.18). 

Finally, we will show that if zo E S^rn H then there exists a solution to (2.17- 
2.18) in the disc DbsL"* (^o)- To this end, we first note that Zo £ J^m H implies that 
(j>{zQ) + 6{zo) is purely imaginary. Combined with the first bound in (5.34) we conclude 
that within a distance of at most (| \og{qm/qn) \ + it)L~'^ + 2e^^^ from cf){zQ), there 
exists a point of the form w = {log{qn/Qm) + i7r(2/c + 1))L~'^ for some integer k. We 
now set e = S2i~'*/2 and 5 = (| logiqm/qn) \+4:)L-^. By the first condition in (5.29), 
we then have |(/)(zo) — w| < S, while the first bound in (5.33) together with the definition 
of B2 implies that S < e|0'(zo)|/2. We therefore can use the first part of Lemma 5.3 
to conclude that there must be a point z' G ^e{zo) such that (t>{z') = w, implying that 
z' is a root of f{z') = 0. Finally, by the already proven statement (3) of the lemma, 
there must be a solution of the equations (2.17-2.18) within a distance strictly less than 
Cie~'^ from z'. Since e + Cie^^ < B2L~'^ by the second condition in (5.29), this gives 
the desired solution of the equations (2.17-2.18) in the disc D^^l-^^C-^o)- n 

Next we will prove Lemma 4.4 which provides a lower bound on f{z) on the bound- 
ary of certain discs. 

Proof of Lemma 4.4. Let a. and M be as in Assumption B3, let k = k/2, and let c 
and -L4 be the constants from Corollary 5.2. We will prove the claim with 



C2 = (2eM||q|U) 



-1 



and C2 = max{c2, 22ea ^} 



(5.35) 



and, given C > C2, with L2 defined by the condition that L2 > L4 and 



CcL < c/L, L'^e 



(I-tL 



< 1, 



(5.36) 
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and 

2e(M + M^)\\q\\^C^L'^eL < 1 (5.37) 

hold whenever L> L2. 

Fix i > Z/2 and choose a point zq e ^k/(4l) n (J^m U J^„) with ^q^^ (zq) C ff. 
Let s < CeL and note that by (5.36) we have s < c/L. Applying Corollary 5.2(1) 
to the disc ^s{zo) we find that DsC-^o) C ^k/{2L) C ^k/l- In particular, the bounds 
of Assumption B are at our disposal whenever z G D^^^(zo). The proof will pro- 
ceed by considering two separate cases depending (roughly) on whether |/(i;o)| is 
"small" or "large." We will first address the latter situations. Let us therefore sup- 
pose that \f{zo)\ > ^L'^eL CC^o)"^ • In this case, we will show that (4.23) holds with 
s{zq) = C2€l. (Note that s{zo) < C^^-l < Ccl by our definition of (72.) A crucial part 
of the proof consists of the derivation of an appropriate estimate on the derivative of /. 
Let s < CeL and let z be such that |z — zo| < s. Recalling the definition (4.2) of bm{z) 
and using Assumptions B2-B3, the second and third bound in (5.36) and the fact that 
one of the values |Cto(^o) | and |Cn(2o) | must be equal to C{zo), we have 

\f'{z)\ = L'^ qmhrnizX^^Xz)"^" +qnhn{z)C\zY 



qmM\Cm{zo)\^ +?„M|C„(^o)|^ 



M\z-zo\L''+L''e-''^ (5.38) 



< 4eM||q|UL'^C(^o)' 



whenever z G ■S^^jh- As argued above, z G D^j^ (-zo) implies that [zq, z\ C so 
by the Fundamental Theorem of Calculus we have 

1/(^)1 > \ j{zQ)\- ^eM\\c^\^L\{z^f\ > 4L'^CW''(eL " ^) (5-39) 

for all z G ]S)s{zo). The bound (4.23) now follows by letting s f C2eL- 

Next we will address the cases with \f{zo)\ < 4L''eLC(^o)^' ■ Let s < CeL and 
pick z such that |z — | = s. This point belongs to the disc ©c^^ (zo) which we recall 
is a subset of S^^jh- The second-order expansion formula 

j{z) = /(zo) + f{z^){z - Zo) + (z - z^f I dt [ dif"{iz + (1 - i)zo) (5.40) 

Jo Jo 

then yields the estimate 

1/(^)1 > |/(^o) + {z- zo)f'{zo)\ - K{CeLfL^''azo)''" (5.41) 

where 

K = Icizoy'-'L-^" sup{\ f"{z)\ : z&^,\z-zo\< CeL}. (5.42) 
Proceeding as in the bound (5.38), we easily get 

K < 2e||q||oo [M'^{1 - L'"^) + ML'"^] , (5.43) 
which implies that K < 2e||q||oo[M2 + M]. 
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It remains to estimate the absolute value on the right-hand side of (5.41). Abbreviat- 
ing bm = bm{zo) and 6„ = 6„(2o), we can write 

f'{zo) = L'^ibmqmC^Kzo)'^'' + br.qnCi^\zof') 

Without loss of generality, let us suppose that |Cm(^o)| > ICn(^o)| and, consequently, 
|Cm(2o) I = ({zo), because zq S =5^™ U Applying Assumption B3 together with the 
assumed upper bound on |/(zo) |, we get 

\{z - zo)f'(zo) + f{zo)\ > (5g™se-^'«"'^ - Asl (1 + sL'^M))L'^C(^o)^', (5.45) 

where we recalled that |z — zo| = s- Since s < Cei, the third inequality in (5.36) gives 
that sL'^M < CML'^tL < 1- Let now s be so large that s>\CeL- Using this bound 
in the first term in (5.45) and using the second inequality in (5.36) we thus get 

|(^-^o)/'(^o)+/(^o)| > {\aC2e-^ -S)L'^eLazot' >^L''^Lazo)''' ■ (5.46) 

Moreover, using the above bound on K and the inequality in (5.37), the last term on 
the right-hand side of (5.41) can be shown not to exceed L'^CLCizo)^ ■ Putting (5.41) 
and (5.46) together with these estimates, we have \ f{z)\ > 2L'^eL({zo)^ for all z E 
^CcL (^o) ^^^^ ^ ~ \z — zo\ satisfies ^Ccl < s < Ccl- The proof is finished by 
taking s t C^l- 

The last statement of the lemma is an immediate consequence of the fact that when- 
ever the above procedure picks 5(2:0) = C2ei, and 62 < C, then the argument (5.38- 
5.39) implies the stronger bound 

inf \f{z)\>2L''eLC{zo)''". (5.47) 

z: \z—zo\<s(zq) 

Now, if / has a root in Dg^g^^ (zq), then this bound shows that we could not have chosen 
s{zq) = C2£L- Therefore, s{zo) must be equal to the other possible value, i.e., we must 
have s{zq) = Ce^. □ 

Proof of Lemma 4.5. We will prove (4.24) with A3 = 2Co||q||i, where Co is as in (2.14) 
for ^ = 0. Let Zq and M be as in Assumption B and let L4 and c be as in Corollary 5.2. 
Let C e (0, 00) and let us choose Z/3 > max{i4, Zo} in such a way that 

max{Ce-^^,Ci''e-5^'T^} < y, MCL^e'''^ < log 2, (5.48) 
+ MCX^^'e-s^''^^ < tL, (5.49) 

and 

1L<^ and MCi^<*e-5^''T^ + iV^-^ < 2d\ogL + log Co (5.50) 
hold for all i > L3. 

We will treat separately the cases zo G "^-y^^ n'^2K/L(^o) and zo € '^^l\'^^2k./l{zo). 
Let us first consider the former case, so that 5l{zo) = e~'^^. The first condition 
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in (5.48), the fact that Dc<5i,(zo)(-2^o) C 6 and 7l < k/{2L) therefore allow us to 
use Corollary 5.2(3), from which we conclude that 'Ocs^(zo){.zq) C '^■y^- For ^ G 
I^C(5i(zo)(-^o) we may thus apply the ^ = version of (2.14) to the function g{z) = 
"{m.?i}.L(^)- Combined with the bound (5.3), the second condition in (5.48) and our 
definition of A3 this immediately gives the desired bound (4.24). 

Next we will attend to the cases when zq Ues in \ so that 5l{zo) = 

2^dg-iLSi.. Letus define Q' as in (5.13) with s = C5l{zq), i.e., 

Q' = {k& n: ^c5,{zo) C y./L{k)]. (5.51) 

By Corollary 5.2(2), the set Q' is non-empty and D(;75^(zo)(^o) C '^k/l{Q')- Let z G 
Dc(5i,(zo)('2o) and let us estimate g{z). We wiU proceed analogously to the preceding 
case; the only difference is that this time we have 

g{z)^SQ,^L{z)+h{z), (5.52) 

where the extra term h{z) is given by 

h{z)= J2 9k[Ci''\z)f. (5.53) 

fc€:Q'\{m,n} 

Now |S'g/.i(z)| is estimated as before: Using that z E ^k/l{Q')^ the bounds (2.14) 
and (5.3) immediately yield that \Sq>^l{z)\ < Co\\q\\iL''^SL{zo)C{zo)^'' . (Here we 
used that the term e'^^'''^^''^''°h~''^ is bounded by e'^^"^^^ < 5l{zo) as follows 
from (5.49).) 

Therefore, we just need to produce an appropriate bound on To that end, we 

note that, since [20, z] C '^k/l{Q') and \z — zo\ < C6l{zo), we have from (5.4) and 
Assumption B2 that 

li'^h^t' < ICf (^o)|'''e^^^''-(^°) < |a(^o)|'''e^^^''^(^°'+^'^"" (5.54) 

whenever k e Q'. Since zq G '^-yl, which implies |Cfe^''(2;)| < C(^o)e~'''^/^ whenever 
k ^ {m, n}, we thus have 

\Ci''\z)f < e^^^'^-^^o^+^'^-^'e-^-^-^'cizof (5.55) 

for every k € Q' \ {m, n}. Using the last bound in (5.50), we conclude that \h{z) \ is 
bounded by Co||q||ii''^L(-2o)C(-2o)^'*. Fromhere (4.24) follows. □ 

5.3. Proof of Lemmas 4.6 and 4.7. Here we will establish the two technical lemmas 
on which the proof of Theorem 2.5 was based. Throughout this section we will assume 
that a multiple point zm G ^ is fixed and that Q = Q(zm). We will also use /, g and ^ 
to denote the functions defined in (4.48^.50). 

Lemma 4.6 is an analogue of Lemma 4.4 from Sect. 4.2 the corresponding proofs 
are also analogous. Namely, the proof of Lemma 4.4 was based on the observation that 
either \f{z)\ was itself large in a neighborhood of zq, or it was small, in which case 
we knew that was large. In Lemma 4.6, the function f{z) is more complicated; 

however, a convenient reformulation in terms of Vandermonde matrices allows us to 
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conclude that at least one among its first {q — 1) derivatives is large. This is enough to 
push the argument through. 

Proof of Lemma 4.6. Abbreviating g = |Q| and using A{q) = 2g39(9+i)/2g!yg and 
the constants K = K{Q) and Lq from Lemma 4.2 and M from Assumption B, let 
e = 1/{3K) and L5 > Lq be such that 

^ML-'Rr. < 2, 2||g||iM« < L'^'^CL and A{q)L'^^-'^/'^eL < e/^/q (5.56) 

for all L > Z/5. A choice of L5 yielding (5.56) is possible in view of (4.51). 

Choosing zq E C, we use F(z) to denote the function F{z) = f{z)£,{zo)~^ . First, 
we claim that if (4.52) fails to hold for some L > L5, then we have 

|i^^^n^o)| < i = 0,...,q-l. (5.57) 

VQ 

Indeed, let us observe that, if (4.52) fails to hold, then there must exist a collection of 
points Zk, with k = 1,. . . ,q, such that 

\zk - zq\ = ^Rl and \F{zk)\ < L'^cl, (5.58) 

for all A; = 1, . . . , g. Further, notice that, for |z — 2;o | < iJi,, we have the bound 

implying \F'-'^Hz)\ < 2J2^^ Q q-m \vm\'^ L'^'^ in view of the first condition in (5.56). In 
particular, we have |F(«)(z)|i?| < 2||q||iM«L-'^ for all z in the i?L -neighborhood 
of Zq. With help of the second condition in (5.56), Taylor's theorem yields 



,9-1 

— — [zk-zo) 



<2L'^€L, k = l,...,q. (5.60) 



Now we will write (5.60) in vector notation and use our previous estimates on Van- 
dermonde matrices to derive (5.57). Let x = {xo,xi, . . . , Xq-i) be the vector with 
components 

g F^^\zo) f Zk - Zo Y . r, 1 , 

xe = Rl J. [-. r , ^ = 0,l,...,g- 1, (5.61) 

and let N = (Nk,i) be the q x g-matrix with elements Nk,i = \zk — zo\^R2^ = {k/qY. 
The bound (5.60) then implies that the vector Nx has each component bounded by 
2L'^eL and so ||Nx|| < i^iL'^e^. On the other hand, since N is a Vandermonde ma- 
trix, the norm of its inverse can be estimated as in (4.8). Namely, using the inequaU- 
ties I detN| > q-iii-'^)/'^ and ||N|| < q, we get 

< < g9(9-l)/2+9(9-l). (5.62) 

But then ||xl| < ||N-i|| ||Nx|| < (i"i^i-^'>l'^2^L'^tL implying 

< £!(L'^i?^)-^||x|| < A(g)L2d-rf/9e^, (5.63) 
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where we used that L'^{L'^Rl) ^ is maximal for £ — 17 — 1, in which case it equals 
jjid-d/q^ -y^i^j^ j^gjp Qf ^j^g YasX condition in (5.56), the claim (5.57) follows for 
all i > L5. 

Having proved (5.57), we will now invoke the properties of Vandermonde matrices 
once again to show that (5.57) contradicts Lemma 4.2. Let y be the g-dimensional 
vector with components 

?An = grne"*'"(^)+"'"(^-^")^'CU)-^', « G Q. (5.64) 

Let O = (0^,m) be the gr x gr matrix with matrix elements Oe^m = v^^- (Here (. takes 
values between and g — 1, while m e Q.) Recalling the definition of F{z), the bound 
(5.57) can be rewritten as | [Oy]^| < e/y^. It therefore implies that 

||Oy|| < e. (5.65) 

The matrix O corresponds to the L ^ oo limit of the matrix M in (4.3) evaluated 
at zm- In particular, since zm G -5^k/l('7i) for all L and all m G Q{zm) and in view 
of the second bound in Assumption B2, the bound (4.5) apphes to O as well. Having 
thus ||0~^ II < K with the constant K from Leimna 4.2, we can conclude that 

||y|| < ||0-i||||Oy|| < K||Oy|| <Ke<^ (5.66) 

using our choice e = 1/{3K). On the other hand, let m be an index for which the 
maximum in the definition of ^{zq) is attained. Then we have 

according to the first condition in (5.56). Moreover, Qm > ^ and thus > ^ in 
contradiction to (5.66). Thus, (4.52) must hold for some s{zo) G [Rl/q, Rl] once L 
exceeds L5. □ 

Lemma 4.7 is also quite similar to the corresponding statement (Lemma 4.5) from 
two-phase coexistence. 

Proof of Lemma 4.7. We will prove the Lemma for Aq = 2e(Co + 3)(M + M^)||q||i, 
where M and Co are the constants from Assumption B. 

Let c and L4 be the constants from Corollary 5.2 for k = k. Since 2;m G ^ is a 
multiple point with Q{zm) — Q, we clearly have that zm G ^e{Q) whenever e is small 
enough. Since i? is open, we also have that {zm) C ^ whenever s is sufficiently 
small. As a consequence, there is a constant Ze = Le{zM) such that zm G ^2k/l(Q) 
'^k/2l{Q) and ^c/Lizu) C & whenever L > Lq. Using Corollary 5.2, we reach the 
conclusion that Ds(zm) C ^k/l{Q) whenever L > max{LQ, L4} and s < c/L. We 
now choose Lq > maxjLg, L4} in such a way that 

p'l<c/L, p'l<2pl, {l + 2pL)e-^^ <{M + M^)pI, ^^^^^ 
4(M + M'^)pIL'^ < 1, e^^"-"^ < 2 

whenever L > Lg. By the above conclusion and the first condition in (5.68), we then 

haveDp'^(2M) C '^k/l(Q) whenever L > Lq. 
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To prove (4.53), let us recall the definitioii of Bq^l{z) in formula (2.13) from As- 
sumption B4. Then we can write g{z) as Eq^l{z)Q{zm)~^ + h{z), where 



meQ 



(5.69) 



Our goal is to show that both l (2^)C(^m) and h{z) satisfy a bound of the type 
(4.53). 

We will begin with the bound on h{z). First we recall the definition of 4>m{L) to 
write 

' c^^ (z) \ / ci^^ (^) \ / d''^ (^m) ^ 



Cm (zm)' 



Cm (^m) 



(5.70) 



The first term on the right-hand side is to the leading order equal to e*"™ ^"'^ , which 
is approximately equal to e""'^^"^")^ . To control the difference between these two 
terms, and to estimate the deviations from the leading order behavior, we combine the 
bound (2.10) with the second-order Taylor formula and (2.11) to show that, for all e 

{zm) and all m G Q, 

|log(C(f)(^)/C(f )(^m)) - vm{z -zm)\< e-^Vi + + M'){p'^f, (5.71) 

where we have chosen the principal branch of the complex logarithm. Combining this 
estimate with the second and third condition in (5.68) and the bound (2.9) from As- 
sumption B2, we get 

\LHog{C^^\z)/C^'^\zM))-Vm{z-ZM)L''-i<f>m{L)\ < ?,{M + M^) pIl" . (5.72) 

Using the fourth condition in (5.68) and the fact that I e"' — 1 1 < e|«;| whenever |tt;| < 1, 
we get 

\h{z)\ < 3e(M + M2)||q||i L'^pl^{zf' . (5.73) 

Now ^{z)^" < ^izo)^''e^'^^^^ < 2^{zo)^^ by the fifth condition in (5.68), so we 
finally have the bound \h{z)\ < A£^{zo)^'^ L"^ p\, with A given by ^ = 6e(M + 

M^)\HW. 

It remains to prove a corresponding bound for Eq^l{z)C{zm)^'" ■ First we recall 
our previous observation that D^'^ (^m) C '^k/l (Q), so we have Assumption B4 at our 
disposal. Then (2.14) yields 



< CoL'*||q||ie- 



-tL 



Ciz) 



((zm) 



zGBp>{zm). (5.74) 



Also, by the definition of '^k/l{Q)^ we have that ({z) = min^gQ |Cm('2)l when- 
ever z e Dp/^(zM)- For z e Dp^(zm), we can therefore find a index m G Q such 
that I Cm (-2) I = C(-2) • With the help of (5 .3) and the bound (2.9) from Assumption B, we 
thus get 



C(^) 



C{zm) 



< 



C{zm) Cm{zo) 



< 



C{zm) 



MRLL-'L-'e- 



(5.75) 
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Combined with the estimate (5.72) for z = zq, and the last three conditions in (5.68), 
this gives 

■ C(£) 1 < ^MR,L<^^L'^e-^^3(M+M-)plL'^^(^^^^L- < 2e^(^o)^'. (5.76) 

Using the third condition in (5.68) one last time, we can bound the right-hand side (5.74) 
by 2eCo||q||i(M + ]\/P)L'^pl^{zo)^\ Combined with the above bound on \h{z)\, this 
finally proves (4.53). □ 
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